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APPENDIX A. TENTATIVE MLS ACCURACY REQUIREMENTS 


GENERAL 

Safe flights and landings generally require the following aircraft per- 
formance : 

o The aircraft path deviations should be within safe limits. 

o The aircraft path rates (such as sink rate) should be within 
safe limits. 

o The aircraft attitude changes should be comfortable to pilot and 
passengers. 

o The aircraft control surface movements should be within reason- 
able mechanical limits and should allow adequate margins for 
response to air turbulence and other factors. 

o The aircraft control column activity should be comfortable to 
the pilot. 

o For manual flight the display activity should be acceptable to the 
pilot. 

The MLS shall not compromise the ability of the aircraft to maintain these 
criteria. 

The factors of error magnitude, duration, spectral content and zone of occur- 
rence are important as well as the factors of aircraft type, AFCS configuration, 
gain, and transient response. 
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Error specifications should be directed toward guidance signal errors which 
are related to these factors, and should register the influence of data rate 
variations and not be unduly affected by irrelevant, higher frequency 
variations. Bias and noise as used by RTCA SC-117 do not account for these 
factors. A method which does is described in this Chapter. 

Path following errors and control motion noise are the concepts used to 
delineate accuracy.* 

DEFINITIONS 

Error and Noise Definitions 

Angular Error -- The angular error is the difference between the processed 
sampled data output and the true position angle at the sampling time. The 
angular error budget is partitioned into two categories, bias and noise. 

Angle Bias (Includes Receiver Bias ) — Bias is the long-term misalignment 
between a specific MLS course and a selected course, and includes the MLS 
ground system as well as the MLS airborne receiver mean errors which cannot 
be reduced to zero by real-time calibration techniques. 


* Since the path following error and control motion concepts for specifying 
accuracy are new, this document also presents the RTCA SC-117 accuracy 
specifications . 
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Angle Noise (Includes Receiver Noise) — Angle noise is defined as spatial 


and temporal perturbations in the guidance signal. It originates from both 
ground and airborne equipment and the environment . 

Path Following Noise -- Path following noise is defined as that portion of 
angle noise which can cause aircraft motion; it exhibits relatively slow 
variations . 

Control Motion Noise — Control motion noise is that portion of angle noise 
which affects control surface, wheel, column motion, and aircraft attitude; 
it exhibits moderately fast variations. 

Extraneous Noise — Extraneous noise is that angle noise which exhibits 
variations too rapid to affect aircraft control and guidance. 

Path Following Error — Path following error is defined as the angular 
deviation from a predetermined course of an aircraft perfectly following MLS 
guidance commands. The error is thus due to angle bias and path following 
noise in the guidance signal. 

Range Error — The range error is the difference between the suitably processed 
DME range and the true distance at any given point in time. DME bias and noise 
errors have the same general definition as for angle guidance. 

Course Linearity Error -- Course linearity error is the deviation of the 
angle coding scale factor from the nominal, about a selected course. Linearity 
errors affect effective AFCS gain and display sensitivity, and con- 
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tribute to aircraft instability. 


Accuracy Zones 

The following zones are defined within the MLS coverage in order to facilitate 
definition of system accuracies. These zones are based on the functional 
use of the MLS data and the operational significance of the MLS errors. The 
zones are as follows: 

o Initiation Zone 
o Maneuvering Zone 
o Landing Zone 
o Roll-Out Zone 
o Missed Approach Zone 

ERROR DETERMINATION 

Path following noise and errors and control motion noise are determined by 
passing the time records through standardized filters. The filter character- 
istics are based on a wide range of existing aircraft response properties, and 
are believed to be adequate for any foreseeable aircraft as well. The fre- 
quency response of the aircraft lateral or vertical/longitudinal channel is 
divided into three major spectral regions--a low, middle, and high frequency 
region, as follows: 

o Low - Aircraft path following components 

(0 to 1.5 radian/second-longitudinal channel) 

(0 to .5 radian/second-lateral channel) 

o Middle - Control surface motions, wheel and column motion, aircraft 
attitude . 

(.5 to 10 radians/second-longitudinal channel) 

(.3 to 10 radians/second-lateral channel) 
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o High - Does not affect aircraft control and guidance. 

In terms of spectral density, the bias would be lower frequency limit of 
noise (approximately 0 to .05 radians/second for a 60-second record). 

While the term "path following error" suggests the difference between a 
desired flight path and the actual flight path taken by an aircraft following 
the guidance, in practice this error is estimated by instructing the test 
pilot to fly a desired course, and measuring the difference between the 
filtered guidance indication and the corresponding position measurement 
determined by a high-accuracy instrument such as a theodolite. The errors 
and spectral distribution thus obtained give an accurate estimate of the 
path following parameters. A similar technique is used to determine the 
control motion noise. 

Treatment of Sudden, Large Errors 

Interference or multipath can occasionally cause large, sudden changes in 
angular indication. Provision shall be made to handle such transients while 
maintaining validation and coast requirements. Capability of "coasting" 
through periods of transients shall be provided, which rejects loss of data for 
a period of time up to 2 seconds, except that those functions, actively in use 
to determine flare altitude, shall be limited to 0.5 seconds coasting time. 

Bias 

The angular bias is determined by averaging the time error record of a test 
flight (the difference between the MLS-derived angle and the tracker-derived 
angle) over a period of 60 seconds. 
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Path Following Error 


Path following error must not exceed the path following error specification 
more than 57. of the time over any 60-second portion of a flight record. 

The flight record here is obtained at the output of the standard path 
following filter. The procedure is described in Figure A-l. 

Control Motion Noise 

Control motion noise must not exceed control motion noise specification more 
than 57, of the time over any 60-second portion of a flight record. The 
procedure is shown in Figure A-l. The flight record here is obtained at 
the output of the standard control motion filter. 

ACCURACY SPECIFICATIONS 

The accuracy requirements for the MLS are presented in Tables A-l and A-2. 

The method of measuring the errors is specified in Figure A-2 and Table A-3. 
The values contained in Tables A-l and A-2 include the effects of multipath 
and EMI. 

The exceptional path following error or control motion noise existing less 
than 5% of any 60-second period which exceed the stated limits shall be of 
such a magnitude and length that they present no hazard to flight or 
excessive strain on the aircraft, its pilot or passengers. 

The azimuth error tolerances are listed in feet at the error window of the 
minimum guidance altitude (MGA) ; the angular error figures are premised on the 
given distances to the error windows. For shorter runways, .the same equipment 
would yield superior guidance (in feet) at the MGA. 
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Figure A-l -- Time Error Record Analysis of Path 
Following Filter Output and Control 
Motion Filter Output 



Table A-l -- Path Following Error Specification, 
Type 3 Equipment 
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Table A-2 — Control Noise Specification 
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Figure A-2 — Allowable Degradation Characteristics 
for Type 3 Subsystem 
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ERROR INTERPRETATION 


Procedures are outlined below which relate the measured parameters to the 

specification for accuracy. Refer to Figure A- 3 for definition of points 

outlined below and to Table A-4 for the filter configurations. 

Point A : MLS Raw Error Data. 

Point B : Time Average over any 60-second Portion of Flight Course. 

Course Bias - See Table A-l for limits. 

Point C ; Path Following Error; see Table A-l for limits; use technique 
described in Figure A-l for calculation. 

Point D : Use Table A-2 for limits. 

Data Rate ; The MLS signal format shall accommodate different data rates 
in different configurations. The format shall be capable of 
providing the minimum information update rates for the functions 
in each configuration as shown in Table A- 3. 

TABLE A- 3. MLS MINIMUM INFORMATION UPDATE RATES 
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Figure A-3 — MLS Measurement Methodology 




Table A-4 — Filter Configurations and Corner 
Frequencies 
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ANGULAR ERROR TOLERANCE RATIONALE 


The MLS error budget uses lineal feet at the MGA as the primary error 
standard. The allowable angular error at the MGA depends on the distance 
of equipment from the MSA. 

Azimuth 

The azimuth angular error figures cited are for the longest runway for which 
that equipment is intended to be used. It should be noted that both A2 
and A3 equipment meet the most stringent requirements in the MGA. The 
degradation factors are structured so as to give highly accurate centerline 
guidance for parallel runway operation (100 feet at 10 NM for type A3 
equipment) , and quite accurate guidance throughout the coverage (600 feet at 
20 NM at ±60° azimuth) for ATC interface and curved path approaches. 

Elevation 

The elevation errors are specified lineally at the MGA. Since the steeper 
slopes intersect the MGA closer to the equipment, the angular error can 
degrade and still maintain the same altitude error (at steeper slopes the DME 
error will predominate anyway). For higher elevation angles at longer range, 
the degradation factors are structured to maintain constant vertical error 
at a given altitude. 
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APPENDIX B. MARSAM II COMPUTER PROGRAM SUMMARY 


INTRODUCTION 

The overall purpose of the MARSAM II program was to develop and implement 
a mathematical computer model for use in the performance assessment of 
reconnaissance sensor systems of varied types operating on prescribed aerial 
flight profiles against ground targets in specified background and weather 
environments. MARSAM II (the mathematical model acronym) addresses those 
aspects of sensor performance as are related to the capability of such 
systems to provide target identification detail. Specifically, the types 
of aerial sensor systems considered in the MARSAM II model are: frame 

and panoramic cameras, television, the visual observer, vertical and forward- 
looking infrared, side-looking and forward-looking radar, and ELINT. As 
applicable to the different sensor types, film record and/or display modes 
of operation are considered. In addition, there is provided as an integral 
part of the MARSAM II computer model a stored library of characteristic data 
for numerous target-elements, backgrounds, and weather conditions (such data 
is readily expanded or modified by the user analyst). Where the provided 
library data base is considered applicable to a given problem, the task of 
preparing model inputs is significantly eased for the user. Available 
outputs from MARSAM II range from detailed sensor system performance para- 
meters and associated probability measures of detection, recognition and 
identification to mission success measures. The MARSAM II computer model 
was developed as a tool for use by sensor systems design analysts in their 
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preliminary sensor performance sensitivity studies and for use by systems 
operations analysts in establishing total reconnaissance or reconnaissance- 
strike system requirements. The usefulness of MARSAM II as an analytical 
tool may be measured in terms of the degree to which the user analyst 
becomes familiar with the extent of analytical detail and assumptions in the 
model and utilizes the model output not necessarily as only an end result 
for his particular problem, but also as meaningful input to additional 
assessment measures he has developed to evaluate system performance. 
References 3 and 11 describe the model and its use in detail. Part 1 
summarizes the objectives, scope, and structure of the MARSAM II mathematical 
model. Part II presents the detailed systems analysis of the reconnaissance 
sensors addressed in the model. 


MARSAM II SUMMARY DESCRIPTION 

Figure B-l illustrates the capabilities implemented in the MARSAM II model for 
the performance assessment of sensor systems of varied types operating on 
prescribed aerial flight profiles against ground targets in specified back- 
ground and weather environments. Specifically, MARSAM II addresses those 
aspects of reconnaissance sensor performance as relate to capability to pro- 
vide target identification detail. In summary, the MARSAM II model has been 
structured for efficient use by sensor system design analysts in their 
preliminary sensor performance sensitivity studies and by systems operation 
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Figure B-l -- MARSAM II Capabilities Summary 
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analysts in establishing or evaluating total reconnaissance or recon- 
naissance-strike system requirements. The particular analytical approaches 
followed and the level of detail considered in the various MARSAM II sensor 
models were dictated to a great extent by the nature and current general 
availability of input characteristic data, particularly that relating 
to sensor performance specifications and to target and background signatures. 
Thus, for example, while (a) the modulation transfer function approach to 
evaluation of reconnaissance sensor system performance and (b) the 
capability to examine in great detail the spectral and spatial variation 
in signatures may be desirable and mathematically tenable, there is a 
lack of sufficient data to exercise such considerations productively. 
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Simplification of model use results from availability of a library of target- 
element*/environment/sensor-system characteristic data, such library 
data base being an integral part of MARSAM. To input a problem to the 
model, the analyst describes target composition by a code number for each 
of the target elements within the target. Similarly, atmospheric and 
background environment characteristics may be input, in the main, by 
specification of code numbers . Thus , such code number inputs are used 
as the means to obtain from the library and input to the model a majority 
of the required target-element, enviroment, and sensor-system characteristic 
data (library data is readily expanded or modified as desired by the 
analyst) . 

Further simplification and efficiency of model use results from the 
automatic manner of sensor-to-target offset consideration. That is, for 
a specified sensor and specified flight speed and altitude, the model 
automatically determines and assesses sensor performance at only those 
aircraf t-to-target offset distances for which at least one target element 
falls within the sensor field of view. 

* In the vernacular of MARSAM, a target is defined as a group of one or 
more target elements where, for example, a target element may be a 
man, a truck, a boat, a hangar, a surface opening, etc. 
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MARS AM PROGRAM STRUCTURE 


As indicated in the summary information flow diagram of Figure B-2, the MARSAM 
program consists of computer models which address the performance of sensors 
generically classed as to operation in the visual, the infrared, or the 
radio-frequency portion of the electromagnetic spectrum. 

The four basic sections of the MARSAM Program are the computer program itself, 
the fixed data deck, the library data deck, and the executive deck. Included 
in the program section are all the required sensor system and subsystem 
option routines. The fixed data section contains data which is constant 
for all sensor models in the MARSAM Program. The library section consists 
of stored characteristic data for target elements, background and weather 
environments, and sensor subsystems. Such library data is available for use 
and modification through code-number call in the execution list. The 
execution list section, generated by the user, selects the data, defines a 
particular type of problem run, and specifies the desired output. 

TV and FLIR Sensor Model Summary 

This section presents a brief description and data flow diagram for the TV 
and FLIR sensor models in the MARSAM II computer programs. The basis for 
development and the detailed analytical treatment for each sensor model is 
contained in Reference 11. Also described here is the Human Factors Display 
Model which is utilized by both the TV and FLIR models. 
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Figure B-2 — MARSAM Program Information Flow 
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Television (TV) Sensor Model (see Figure B- 3) •" 

The television sensor model addresses near real-time forward -viewing airborne 
television sensors. Alternate modes of system operation addressed within the 
model include (a) moving display, (b) stationary display, or (c) tracking. 

The moving display mode refers to a continuous and real-time presentation 
of data to an observer. The stationary display mode refers to the selection 
of a single frame of data and display of that frame for a period of time 
which equals or exceeds the normal frame time. While the two previously 
described modes generally apply to fixed sensor depression angles, the 
third mode refers to target-element tracking by continuously changing the 
sensor viewing angle so that the element is maintained in the field-of-view 
center until overflown. The model attempts only to predict the observer's 
target-element recognition performance under conditions of tracking; i.e., 
tracking accuracy is not evaluated. Outputs include probabilities of 
target-element detection and recognition. 
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Forward-Looking Infrared (FLIR) Model (see Figure B-4) 


The forward-looking infrared model is addressed to forward-looking scanning" 
systems in which imagery from successive scans is presented in near-real 
time to an operator-observed display. Alternate modes of operation considered 
in the model include (a) moving display, (b) stationary display, or (c) track- 
ing. System operation within the 0.3- to 15-micron spectral region is 
considered. Model outputs include predictions of acquired target-element 
detection and recognition probabilities. 

Human Factors Display Model (see Figure B-5) 

The human factors display viewing model is drawn from model developments 
performed by Defense Research Corporation (DRC) under the sponsorship of 
the Advanced Research Projects Agency (ARPA) . Depending on the sensor 
system considered, up to five probability submodels are used to determine 
target-element detection and recognition probabilities, i.e., an element 
size and contrast submodel, a search mode submodel, a confusing objects 
submodel, a signal-to-noise submodel, and a recognition submodel. 


B-10 


11 II I I M II I 1 H I 1 II I I 1 £ 





3 

o 

I- 1 

co 

4-1 

CO 

Q 


QJ 

XJ 

O 




TJ 

0) 

t-r 

CO 

U 

d 


00 

c 

o 

o 

hJ 

X) 

CO 

g 

o 


<r 

i 

« 


<v 

u 

Z> 

00 


B-ll 


u i a a a u i m m i a a a a a a a a a 











it i i i i m i m i a i i i a i i i £ 


Figure B-5 — Human Factors Display Model Data Flow 




LIBRARY DATA 


This section contains and references the characteristic environment, 
target-element, and sensor subsystem data stored in the MARSAM library. 

The ten major categories of stored data, accompanied by a brief description 
of the types of data contained, are listed below: 
o Terrain Characteristics 

Line-of-sight probability data for six terrain types 
o Weather Characteristics 

Data applicable to the photographic, television, visual 
observer, and infrared sensor models for five weather 
conditions 

- Data applicable to the radar sensor models for ten weather 
conditions 

o Turbulence Characteristics 

- Data applicable to the photographic sensor models for three 
turbulence types 

o Target-Element Signatures 

- Dimensions, photo/visual reflectivity, emissivity, and 
radar cross-section for 81 target elements 

o Background Signatures 

Photo/visual reflectivity, emissivity, and normalized radar 
cross-section for 15 backgrounds. 

o Sensor Subsystem Characteristics 

Performance characteristics for: 

A forward-looking infrared subsystem 

A TV image orthicon subsystem 
A TV vidicon subsystem 
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A forward-looking radar subsystem 
A side-looking radar subsystem 
A panoramic camera subsystem 


o Lens Subsystem Characteristics 

Performance characteristics for two lenses, one applicable 
to the photographic sensor model and one applicable to 
the television sensor model. 

o Film Subsystem Characteristics 

Data applicable to photographic systems for three film 
types and three developer/development-time combinations 
for each film 

Data applicable to vertical infrared systems for one film 
type with seven developer /development- time combinations 

o Filter Subsystem Characteristics 

- Filter factor and filter function data for four filter 
types 

o Display Subsystem Characteristics 

- Performance characteristics for one cathode ray tube 


Computer listings of all library data are given in Reference 3. 
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APPENDIX C 


SCATTERING OF ELECTROMAGNETIC WAVES 

INTRODUCTION 

Various sensors which could be used for an ILM depend on their operation 
on the reflective properties of the terrain in the microwave region. 

Many measurements of these properties have been made, however, the pre- 
ponderance of the data is at incidence angles between 10 and 80 degrees . 

Since ILM sensors must operate with incidence angles from about 84 to 89 
degrees, these data are not directly usable. 

In an attempt to obtain usable data, a theoretical formulation of electro- 
magnetic scattering was developed and programmed on a computer under the 
assumption that if the theoretical model corresponded to the measured data 
at those points where data was available , the model could be used to generate 
the needed data at higher incidence angles. The model used was a statistically 
rough surface using the Kirchoff approximation and at first appeared to give 
good correspondence. It was later noted that the formulation was 
missing a cosine of the incidence angle. After correcting the model, no 
set of parameters in the theoretical model could give correspondence with the 
measured data. Several explanations are possible for this. It has been pointed 
out that the source from which the model was obtained has an error in the 
dominant term of the equation at high incidence angles. Another problem is that 
only diffuse reflections are considered in the model with specular components 


C-l 


U H I I I H I I 1 I If I M H H I I I I 



handled independently. When the incidence angle is very high, the diffusely 
backscattered radiation may be less than the specularly backscattered 
radiation. Specular backscatter at high incidence angles seems a contra- 
diction in terms but for a finite illumination it is meaningful. Since the 
illuminated surface is determined by the antenna pattern, this component 
depends primarily on the antenna used to measure it. 

Because of the lack of success in obtaining a good theoretical model for 
scattering, the model, with the cosine removed, was used as an empirical 
fit in some sensor analyses. The total exercise demonstrate the need 
for subsequent measurement of typical terrain scattering at these angles. 

SCATTERING PRINCIPLES 

When an electromagnetic wave strikes an object, it is partially reflected 
in all directions. This reflection is the basis of radar, and extensive 
measurements have been made of the reflection from natural and man-made 
surfaces for use in radar-performance calculation. However, most of 
these measurements have been made in the backscattered direction, at 
incidence angles between 10 and 80 degrees. Therefore even for radar sensor 
analysis at high incidence angles the data is not directly usable. 

The scattering in other directions is also important for sensor analysis. 

It is the source of multipath distortion, and a major contributor to mutual 
interference between sensors operating from different platforms. To 
properly evaluate sensor performance, it is therefore necessary to establish 
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the scattering properties of the surfaces involved. 


The technique selected to establish the scattering properties is to obtain 
a theoretical model of scattering which can be justified from physical 
principles. Any such model will have various parameters which can be 
adjusted to modify the scattering properties, so that it can represent 
various surfaces. The model will be applied to the measured backscattered 
data to establish the values of the parameters for the surface, and these 
parameters will then be used to generate angular relationships. 

Definitions for Scattering Parameters 

The most physical scattering parameter is the reflection coefficient. 

It is defined in terms of the field strengths at the reflecting interface 
as 



where R is the reflection coefficient 

Ej^ is the reflected field 

Ej is the incident field 

While this is an adequate parameter for perfectly conducting, finite 
size, smooth targets, it is difficult to use. 
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A more tractable scattering parameter is the differential radar cross 
section, derived from the radar range equation. The radar range equation 
is 


P t Ct 


4 7i r'' 


4 TT 
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2 

A G r 
4 TT 


(C-2) 


where 


P r 

is 

G t 

is 

a 

is 

A 

is 

% 

is 

r 

is 


the received power 

the gain of the transmitting antenna 

the radar cross section 

the wavelength 

the receiver antenna gain 

the range from the radar to the target 


The first term expresses the power density at the target. The second 
term is an expression which relates scattering to an equivalent area which 
captures all the incident power and isotropically reradiates all the 
power. The third term represents the capture of all the energy impinging 
on the receiver aperture . 


Since the scattering of an isolated target is neither perfect nor isotropic, 
the radar cross section is a complicated function of incident and reflected 
angles, frequency, polarization, target geometry, etc., and contains all 
the deviations from perfect isotropic scattering. 


C-4 


u i i r a a i i i i i i i a a a k i i 



For an area target, the concept of cross section is generalized to a 
dimensionless cross section per unit area, such that the radar range 
equation becomes: 


P t A 


//( ~3 * 
(41 Y) r 


y G t G r dS 


(C-3) 


where s is the reflecting surface 


Y is the differential cross section per unit 
surface area. 

G t ,G R are the antenna gains in the direction of the 
differential area dS 


This equation is often used in the form: 

p - P t G t e R > 2 fr ds (c _ 4) 

1 < 4„> 3 r 4 > 

which assumes that the gain across the antenna beamwidth is constant 
and zero outside the beam. By examining this equation, it can be seen 
that for finite targets 

y = r 2 (C -5) 

A 

where A is the target area 

R is the reflection coefficient 

A is the wavelength. 
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The scattering geometry which will be used throughout is shown in 

Figure C-l. Horizontal polarization is defined as having the E vector in the 

XY plane, vertical polarization has the E vector in the XZ plane. 


SMOOTH SURFACE REFLECTION 

Scattering from an infinite smooth surface is the familiar mirror 
reflection obeying Snell's law. That is, any energy that is not absorbed 
or transmitted is reflected such that the angle of incidence is equal 
to the angle of reflection. Furthermore, the amplitude of the reflected 
field strength and the change in phase at reflection depends 


C-6 


ii I I i I I 1 H I I If I I I li Jf fi l £ 




Figure C-l — Scattering Geometry 
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on the Fresnel reflection coefficient, which is a function of the 
polarization of the incident wave and the properties of the surface. 
For non-magnetic surfaces, the coefficients are: 


K « 


: cos 0. -A# e - sin^ 0. 

_r 1 V r i 

: cos 0. +"\/e - sin^ 0. 

r i y r i 
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cos 0 . + a e - sin^ 0. 
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where 0^ is the angle of incidence 

e r is the complex permittivity 

r v is the reflection coefficient for vertical 
polarization 

R H is the reflection coefficient for horizontal 
polarization 


For a perfectly conducting surface, e r becomes infinite and Rv=1=-Rjj. 

Thus, all the energy is reflected, the vertically polarized portion with 
no phase change and the horizontally polarized component with a 180° phase 
change. However, if the surface is not perfectly conducting the reflection 
coefficient at low grazing (high incidence) angles becomes Rv=Rh=- 1. Some 
values of the reflection coefficients are shown in Figure C-2 for typical 
earth and sea permittivities. 
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If the size of the surface is less than infinite, the energy will not be 
redirected perfectly due to the effect of fields on the edges of the surface. 
As the surface becomes smaller, the spread of the reflected beam increases. 
This aperture phenomenon is the szme as the diffraction in optics, or the 
gain functions for reflecting or aperture antennas. Figure C-3 illustrates 
this phenomenon. 

The differential cross section of the target is: 

2 2 2 
Ypg = XY sine tt_Cx- X_ sine tt Ky Y |ct pg | ( C _ 8 ) 

X XX 

where X is the dimension in the x direction 

Y is the dimension in the Y direction 
sinx 

sincx = “ 

= sin 0i - sin 02 cos 0 
=-sin 0 2 sin 0 

a is the polarization dependent reflection 
coefficient 

g is the transmitting polarization 
p is the receiving polarization 
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For linear polarizations. 


a yy = -cos 02 Cos 0 Ry 

(C-9) 

a m = sin0R v 

(C-10) 

a vj[ = cos ©jcos 0 2 sin0Rn 

(G-ll) 

a HH = cos ®1 cos0Rjj 

(C-12) 


STATISTICAL SURFACE MODELS 

Most scattering of interest does not involve smooth flat surfaces, but 
rather rough surfaces. Grass, concrete, weeds, lake or ocean water are 
all examples of rough surfaces. A high sea state is considerably rougher 
than a concrete slab, but they are both rough. 

Equations 

The equations for scattering can be derived by assuming that the electric 
field at any point on the surface is the field that would be present on 
a tangent plane to that point. Then, either the Helmholtz integral can 
be evaluated in the far field (Ref. C-l) or the total area of "specular 
points" which have parallel tangent planes can be evaluated from the 
statistics of the surface roughness (Ref. C-3). Either technique leads 
to the same formula for scattering. 

The Helmholtz integral at a point in the far field distance r from the 
surface is (Ref. C-l) 

E - j exp (jKr) f iv r 

4 tt ar Js (EV-p) *n e * dS (C-13) 

where E is the electric field strength at the point pe^ 
unit incident field 

F is the Fresnel reflection coefficient 

V = k 1 -k 2 

C— 12 


li a i a a a a a a a a a a a a a a i a 



K]^ is the propagation vector of the incident wave 
K2 is the propagation vector of the reflected wave 
n is the unit normal to the surface 
p = + K2 


For a finite surface area in cartesian coordinates, 
jk exp(jkr) 


E = 


4 11 r 




+c -b) 
x y e 


> > 

jv*r 


dS 
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where a = ( 1 -F) sin© ^ + ( 1 +F) sin ©2 cos 0 

b = ( 1 -F) cos 0 2 - ( 1 -F) cos 

c = ( 1 + F) sin ©2 cos# 

p (x,y) is the height of the surface at (x,y) 

X dx 

Since the field due to a perfectly conducting area of the same size is 


jkA Cos ©^ exp (jkr) (C- 15 ) 

2 7T r 

A reflection coefficient for the rough surface can be defined as 

1 C .>.> 

p = Js(a#' +C q' -b) e lv * r dS (C- 16 ) 

2ACOS©! X y 

where > 2tt r, „ > .> >■_, 

v = — l_(sin 0^-sin©2 cos 0 ) i-sin©2sin0j-(cos©^+cos02) K J 

i,j,1c is the unit vector triad 


C -13 


u a 1 i u u 1 u u 1 11 a i 11 it x a 1 1 



F is the local reflection coefficient. Neglecting the effects of fields 
at the edges of the surface and assuming the surface to be isotropically 
rough, the mean value for a perfectly conducting surface is (ref. C-l ) 

<P> = X(V Z ) sine V x X sincV Y (C-17) 

“Z - T~ 

where V x is the x component of v 

X (v x ) is the characteristic function of the random variable § 


The mean value of pis thus the reflection reduction factor to account 
for the finite size of the reflecting surface. Its value is of significance 
only near the specular direction (incidence equal to reflection) , where 
the definition of near depends on the size of the surface. 


The variance of P has values of significance at angles other than specular. 

O 

Since the received power is proportional to p , power will be scattered 
over other angles. This power is referred to as the diffuse reflection. 

The mean square value of the reflection coefficient, obtained by multiplying 
p by its conjugate, taking the expectation, changing to cylindrical 
coordinates and integrating overall angles is (Ref. C-l ). 
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< e tV Z TdT 
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where F is the local Fresnel reflection coefficient defined on 

pp. 020 and 021 
V 

xy 


4 


2 2 
V + V 
X y 


J o is the Bessel function of order zero 
T is a dummy radial distance variable 
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which 


The integral in the expression for <pp*>depends on < e^ v t ^ 

is the characteristic function of the probability distribution of % . 
Therefore, to proceed further it is necessary to define this distribution. 

By the central limit theorem, it can be expected that the surface is 
Guassianaly distributed in height. Two correlation functions, corresponding 
to "peaky" and more smoothly bumpy surfaces will be investigated. The 


characteristic function for a gaussian surface is: 

2 

< e iVz( -S'-5) > = ex P (-v2 z a (i - C (t) ) 


(C-19) 


where a is the variance of surface height 
C(t) is the correlation coefficient 


.2 2 


Slightly Rough Surfaces 

If the surface is slightly rough in the sense that a 1 " < < 1 t 
then the power series for the exponential will converge rapidly enough that 
only the first term has significant contribution to reflection. In this 


case, the integral is readily evaluated yielding (Ref. C-l , C-4) 


ff F 
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_ „ 2 2 Vxy T 
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if C (t) =e T 


2 2 2 2 
< 2tt|f| Vz g T exp (-Vz 2 q 2 ) 
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The first form corresponds to smooth bumps, the second form to peaks. In 
both cases the specular component is ignored. 


The differential cross section, obtained by multiplying by the target 
aperture gain and rearranging terms is 

.2 


y =(2 ir ) 4 (|) 2 (y) 2 | F | 2 (^) 2 Exp (2ir Y 
XJ.2 , V z.2 , ,T.2 ,V . I! 

<T> + <x> J 


for C(r) = e 


2.2 

- x /T 


O/O x4| | 2 / V z, 2. a .2/,1 .2 _ ,9 N 2, V Zv2 

= 2(2ir) I F | (— ) ( — (— ) exp-27r(— ) (— ) 


a + ( 2 ,> 2 (£ 2 ( - f ) 2 > 3/2 
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(C-23) 


for C (x) = e ' T 


/T 


In this form the independence of the basic equations from frequency is 
evident. Except for defining roughness and correlation length in terms 
of wavelengths, these equations are independent of frequency. 
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Very Rough Surfaces 


The other case which can be easily evaluated is the very rough surface. 


in the sense that V a > > 

z 


2 

1. Since v a is very large, the characteristic 
z 


function will have significant values only near t= 0. Further, since the 
correlation function is by definition an even function, it's McLaurin 
series will contain only even powers. Since Cj(0)=l, this term contributes 
nothing to the expansion of the characteristic function about zero. 

Thus the only term of interest is the second derivative. Taking this 
expansion and performing the integration (ref C-4). 
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Converting to differential cross section. 
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for C(t) = e l T l /T 
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where the equations are again written to emphasize the frequency independence 
of the basic equations. 

Shadowing 

At very low grazing angles, significant portions of the surface may be 
shadowed by the surface roughness, as the back sides of hills are shadowed 
from the sun. Obviously, these areas which are in the shadow cannot 
contribute to the reflected power. 


To a first approximation, the effect of shadowing can be considered by 
multiplying the differential cross section by a shadowing function which 
represents the probability that a ray path is blocked from hitting an 
area of the reflecting surface. For backscattering from a normally 
distributed surface, regardless of surface correlation properties, this 
shadowing function is given by (ref C-5) . 

Sj = exp. (-1/4 tan ©x erfc(Rcot ©i)) (C-28) 

where ©x is the incidence angle 
T 

R = for gaussian correlation 

for exponential correlation. 

\f2a 

The function S is shown for various values of k in Figure C-4. If 
kcot 9^>>1 , the complementary error function can be approximated by its 
asymptotic expansion yielding: 


Sx = exp (- 




2 2 

2 n -R cot 0 . 
tan 9 e 1) 


(C-29) 
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Figure C-4 — The Shadowing Function 



For bistatic scattering, the only change necessary is to allow for the 
probability that an area which is illuminated cannot be seen by the receiver. 
If the reflected angle is less than the incident angle, i.e., nearer the 
zenith, this probability is zero. Otherwise, the shadowing function is 
of the same form as the incidence shadowing function, replacing the 
incidence angle with the reflected angle. Thus 
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(-1/4 tan 02 er f c ( Rcot 02^ 
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The overall shadowing function is the product of the incidence and 
reflected shadowing functions or 

3 = S I S R (C-31) 


Polarization Dependency 

The local Fresnel reflection coefficient used in the scattering equations 
is a function of the angles of incidence and reflection, the local slope 
of the surface, the complex permittivity of the surface, and the polarization 
of the incidence and received waves . 


Assuming that all reflected energy comes from specular points, the local 
Fresnel coefficient for arbitrary scattering angles and linear polarization 
can be derived by geometry (Reference C-7): 


F W = 


a 2 a 3 R v (i) + sin02sin 0R H (i) 
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(C-32) 


C-20 


K 


U U MM MM MEM I I I I II 1 i 1 



(C-33) 



(C-34) 

(C-35) 


where = l+ag 

&2 = cos 6].sin ©2+sin 01 cos 02 cos0 
ag = sin 0-^ cos 02+cos 0^ sin 02cos0 
*4 = cos 0^ +cos ©2 
a^ = sin O^sin 02cos0-cos O^cos ©2 
i = arc tan 

is the local Fresnel coefficient for J polarized 
transmission and K polarized reception. 

Since this set of equations defines a scattering matrix, the reflection 
for arbitrary polarizations can be derived from it. In particular, for the 
circular polarizations. 
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For linearly polarized backscattering , the expression for F w and Fgjj 
become the same since i goes to zero. Thus, for backscatter, 

Fvv - ^ V°> ' R * M (C ’ 40 > 

a l a 4 cos0. 

1 

Since measurements have shown differences between horizontally and vertically 
polarized backscatter, some other factor must cause the difference. A 
plausable explanation lies in the assumption of isotropic roughness. Consider 
the bump of Figure C-5. If radiation is incident on this bump and 
specular point backscatter occurs, the field at the surface appears to 
be locally horizontally polarized, regardless of its polarization relative 
to the mean surface plane. However, the slope or correlation distance 
seen by a horizontally or vertically polarized wave can be considerably 
different. In the figure, a horizontally polarized wave sees a steeper 
slope than a vertically polarized wave. Therefore, it is possible that 
different surface parameters should be used for different polarizations. 

Computer Program 

Since the scattering equations of the previous sections are extremely 
complex, a FORTRAN subroutine has been written to compute the differential 
cross section. All of the models given in the previous sections are 
included in this program, as are all combinations of vertical, horizontal, 
and circular polarization. 
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Figure C-5 — Isotrophic Roughness 



Being a subroutine, the scattering equations can be included in other analysis 

programs. The subroutine is called by the statement 

CALL STSCAT (AINC, AREFL, AOAZ, PRMTIV, IPOL, ICOR, ROUGH, SLOPE, GAMDB) 

where AINC is the incidence angle ©i, in degrees 

AREFL is the reflection angle, ©2, in degrees 

AOAZ is the angle off azimuth, 0, in degrees 

PRMTIV is the surface permittivity, a complex number 

IPOL is an indicator of polarization, 

IPOL = 1 vertically transmission and reception 

IPOL = 2 vertical transmission, horizontal reception 

IPOL = 3 horizontal transmission, vertical reception 

IPOL = 4 horizontal transmission and reception 

IPOL = 5 left circular transmission and reception 

IPOL = 6 right circular transmission and reception 

IPOL = 7 right circular transmission, 
left circular reception 

IPOL = 8 left circular transmission, 
right circular reception 

ICOR is an indicator of surface correlation type 
ICOR = 0 gaussianly correlated surface 
ICOR = 1 exponentially correlated surface 
ROUGH is the scale surface roughness a/ 

SLOPE is the scale surface slope 1/ 

GAMDB is the differential cross section y, in decibels. 
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A main program, used to test the subroutine and establish surface 
parameter values has also been written. The test program obtains back- 
scattering cross sections from 10 to 80 degrees, and plots these cross 
sections along with a set of values which are input. By examining the 
resultant plot, it can be immediately determined how well the theoretical 
results match the measured data. 

The name of the test program is TSTSCAT. The program reads in from the 
teletype values for the constant (measured) plot, scale roughness, scale 
slope, and permittivity. Polarization and correlation are changed by 
modifying the program. The call to STSCAT is made with continuation lines 
to enable modifying IPOL and ICOR with literals. Line 180 is IPOL, line 
190 is ICOR. Thus, to enter a value of IPOL, the following program 
modification would be made 
180 &[lP0L] , 

where ^IPOLj is the value of IPOL, to change correlation, the 
equivalent modification is 
190 & [iCOR] , 

A listing of the entire program, including the test drive program and 
the subroutine STSCAT follows: 
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1**RUN »;LlBRARY/PLOT ,R 

010 DIMENSION RESULT ( 10) « ANGLE (10) » I NGAM ( 10 ) , YPLT ( 2 > 

011 REAL INGAM 

12 REAL KPLT 

020 CHARACTER REP*3«YES»3 

021 COMPLEX PRMTIV 

23 YES= "YES" 

24 REP="N0P" 

030 PRINT 60 

031 31 PRINT 34 

032 READ 33. (lNGAM(l) .1=1.10) 

033 33 FOPMAT ( 1 0F6 • 1 ) 

034 34 FORMATI’INPUT ARRAY TO MATCH.") 

040 40 PRINT 70 

050 READ 80.RGH 

060 60 FORMAT("INPUT PARAMETERS.*) 

070 70 FORMAT (-SURFACE ROUGHNESS") ' 

080 80 FORMAT ( F6 • 4 ) 

090 90 FORMAT ("SURFACE SLOPE") 

100 PRINT 90 

110 READ 80.SLP 

111 PRINT 112 

112 112F0RMAT ("PERMITTIVITY".) 

113 READ 114, PRMTIV 

114 114F0RMAT (2F6.1) 

115 PRINT 115, PGH 

116 l 1 5F0RMAT ( / / "FOUGHNESS= " , F6 . 4) 

117 PRINT1 18, PRMTIV 

118 1 18F0RMAT ( "PERM I T I V I T Y r " , F6 . 2 , " ♦ J" ,F6.2) 

120 DO 10 1=1,8 

130 THET A1 = 1/ 5.729 

140 THETA2=-THETA1 

150 PH 1=0 

160 CALL STSCAT(ThETA1,THETA2«PHI, 

1706 PRMTIV, 

1806 4, 

1906 1, 

2006 RGH. 

2106 SLP, 

2206 SIG) 
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230 RESULT < I ) = S I G 

240 10 ANGLE (11=10*1 
2 SO PRINT 260. SLP 
260 260F0RMAT ( " SLOPE = " «F6 • 2 ) 

270 PRINT 20. (ANGLE(I) .1 = 1.10) 

280 PRINT 30. (RESULT ( I ). 1 = 1.10) 

281 PRINT 282 

282 282FORMAT (//) 

290 15 CONTINUE 
300 50 CONTINUE 

305 PRINT 306 

306 306FORMAT ("PLOT ") 

307 READ 380. REP 

308 IF(REP.NE.VES) GO TO 338 
310 20 FORMAT ( 1 OH INC I DENCE .10F6.2) 

320 30 FORMAT ( 1 OHCROSS SECT.10F6.2) 

321 NMPT=2 

322 YMAX=- 10 

323 YMIN=-50 

324 CALL PLOTIXPLT.YPLT.YMAX.YMIN.NMPT.1.36) 

325 XPLT=0 

326 LPLT=1 

327 DO 337 JPLT=1.36 

328 XPLT=8*2»JPLT 

329 LPLTl=LPLT.l 

330 YPLT (1) =KPLT/5.0* (RESULT (LPLT1 ) -RESULT (LPLT ) ) .RESULT (LPLT, 

331 YPLT (2 >=KPLT/5.0»( INGAM (LPLT1) -INGAM (LPLT) ) ♦ I NGAM ( LPLT ) 

332 CALL PLOTIXPLT.YPLT.YMAX.YMIN.NMPT. 0.36) 

333 KPLT=KPLT*1 

334 IF(KPLT.LT.5)G0 TO 337 

335 KPLT =KPLT-5 

336 LPLT=LPLT*1 

337 337CONTINUE 

338 338PRINT 370 

340 READ 380. REP 

350 IF (PEP.EO.YES)GO TO 40 

361 PRINT 364 

362 READ 380. REP 

363 IF (REP.EO.YES)GO TO 31 ‘ 

364 364F0RMAT ("ANOTHER ARRAY".) 

365 370 FORMAT ("ANOTHER CASE") 

370 380F0RMAT (A3) 

380 STOP 

390 END 

1010 SUBROUTINE ST SC AT (AINC .AREFL . AOAZ .PRMTI V. JPOL . I COR .ROUGH. SLOPE .GAMDB) 
1020 COMPLEX PRMTIV.BVV.BVH.BHV.BHH.RVERT.RHORZ.BETA.BCONj.CTFMP 

1030 SM I NC = S I N ( A I HC ) 

1040 SNREF=SIN(AREFL) 

1050 SN0A2=SIN!A0AZ) 

1060 CSI.NC = COS(AINC) 

1070 CSREF=COS (AREFL) 

1080 CSOAZrCOS(AOAZ) 

1082 TEMP1=AINC*AREFL 

1083 IF(TEMP1.LT. .00001) GO TO 430 

1090 A5=CSINC»CSREF-SN!NC»SNREF*CS0AZ 
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1100 TEKP1=(1-A5)/U»A5) 

1110 TE.MP1 = SGRT (TEMPI) 

1120 L ! NC: AT AN ( T EMP 1 ) 

1130 SL INC: S IN (L INC) 

1 1 AO CL I NC = COS (L INC ) 

1150 CTEMPrPPMT I V-SL INC**2 

1160 CTEMPrCSfiRT (CTEM?) 

1170 RVERT: (PRMTIV*CL INC-CTEMP) / < PRMT I V*CL I NC«CTEMP) 

1180 RHORZ: (CLINC-CTEMP) / (CLINC*CT£MP) 

1190 A2 :CS I NC* SNPEF ♦ SNINC*CSREF*CSOAZ 

1200 A3:SHINC»C.SREF*CSINC»SNREF*CS0AZ 

1210 IFUPCL.GT.AIGO TO 230 

1220 GO T0(230.270.310,350) ,IPCL 

1230 230 BVV:A2*A3*RHORZ*5NI.NC*SNR£F»SNOAZ**2«RVERT 

1 2 A 0 1F(1PCL.GT.A)G0 TO 270 

1250 BET A = BVV/ ( 1-A5 ) 

1260 GO TO AAO 

1270 270 BVH:SNOA7*<Sf':REF*A2»»VERT-SNINC*A3*RHORZ) 

1280 IF ( IPCL.G1 . A) GO TO 310 

1290 , BETA:BVH/(1-A5) 

1300 GO TO AAO 

1310 310 eHV:SNOAZ*(SNREF*A2»RHORZ-SNINC*A3*RVERT) 

1320 I F ( I PCL «GT • A ) GO TO 350 

1330 BETA=BhV/ ( 1-A5) 

l 3A0 GO TO AAO 

1350 350 BHH:- (SNINC*SNREF*SNOAZ*»2*RHORZ*A2*A3*RVERT) 

1360 IFUPCL.GT.AIGO TO 390 

1370 BETA=BHH/ (1-A5) 

1380 GO TO AAO 

1390 390 I F ( I POL • LE.6) GO TO A20 

1 AOO BETA: ( (REAL (BHH) *REAL ( BVV) -A I MAG (BHV) ♦ AIMAG (BVH) )**2* (ATMAG(BHH) 

<A 1 0 6 -RFAL (BVH) ♦AIMAG (BVV) ♦dfai. (BmV) > **2 > /A 

IA20 GO id **50 

1A30 A20 BETA: < (REAL (8HH) -REAL (BVV) -AIMAG (BVH) -AIMAG (BHV) )**2* (AIMAG (BHH) 
1AA00 ♦PEAL (BVH) -AIMAG (BVV) »REAL (BHV) ) **2 ) /A 

1 A50 GO TO A50 

1 A58 A30 BETA=(PRMTIV*1-2*CS0RT(PRMTIV))/(PRMTIV-1) 

1 A60 AAO BETA:BETA«CONJG(BETA) 

1A70 A50 EXYSQ=SNINC»*2*SNR£F**2-SNINC*SNREF*CSOAZ 

1A80 IFUCOR.GT.OIGO TO 580 

1 A90 I F ( (ROUGH* (CS I NC*CSREF ) )**2»LT.l) GO TO 5A0 

1500 , 6AMMA:2»BETA* (SLOPE/ (ROUGH* <CSINC»CSREF>*«2) ) **2*EXP (-.25* (SLOPE 

15106 / (R0UGH*(CSINC«C5REF) ) ) »*2«EXYS0) 

1520 X = SLOPE/ (2*R0UGH) 

1530 GO TO 670 

15 AO 5 AO GAMMA: 1.5585A55E3*BETA*( SLOPE* ROUGH) **2*EXP(-39.A78a17* ( (ROUGH* 
1550 6 (CSINC*CSREF) )**2*.25*SL0PE**2*EXY50) ) 

1560 560 K:SL0PE/(2*R0U6H> 

1570 GO TO 670 

1580 580 IF( (ROUGH* fC5INC*CSREFJ 1**2. LT. I ) GO TO 630 

1590 TEMPI :39.A78A17*( (CSI NOCSREF) ‘ROUGH ) **A/SLOPE««2*EXYSO 

1600 GAMMA: 1 2 . 56637*BET A*ROUGH** 2 / ( SLOPE* TEMP l* SORT (TEMP l ) ) 

1610 K = SLOPE / ( 1 • A1 A*ROUGH) 

1620 GO TO 670 

1630 630 TEMPI = 1*39. A7A817*EXYSO*SLOPE**2 
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1640 
1650 6 
1660 

1670 670 

1680 

1690 

1700 

1710 

1720 

1730 730 
17*0 

1750 750 
1760 760 
1770 
1780 
1790 
1800 

1810 810 
1820 820 
1830 830 
18*0 

1850 850 
1860 

1870 870 
1880 880 
1890 890 
1900 


GAMMA- 3 ,1170909E3*8ETA* (ROUGH ‘SLOPE )**2*EXP(-39.*78*17»{ ROUGH* C 
CS I NC*C5REF ) ) *«2> / (TEMP 1*SQRT (TEMPI) ) 

K : SLOPE/ ( 1 • * 1 **POUGH) 

IF (CSINC.GT.CSPEF)GO TO 750 
TEMP1=A0S(SNREF/C5REF) 

IF(K/TEMP1.LT.1)G0 TO 730 
TEMP2 -EXP (- (K /TEMP 1)**2) 

SHADR = EXP (-TEMP1**2*TEMP2/(3.5 5*9*IC) ) 

GO TO 760 

SHADR=EXP<- (TEMPI -K) /3. 55*9) 

GO TO 760 
SH ADR - 1 

TEMP 1=ABS<SN1NC/CS INC) 

IF(K/TEMP1.LT.1)G0 TO 810 
TEMP2=EXP (- (R /TEMPI )**2) 

SHAD I =EXP<-TEMP1*«2*TEMP2/ ( 3 . 55990707*K ) ) 

GO TO 820 

SHAOI=EXP<- (TEMP 1-K 1/3.5**9077) 

GAMMA = GAMM A * SH ADR* SHAD I 
I F (GAMM A. GT « 0) GO TO 870 
PRINT 8 80. GAMMA 
GAMDB-99.99 
GO TO 890 

GAMDBslO»ALOG 10 (GAMMA) 

FORMAT ( 11HGAMMA ERROR.F6.2) 

RETURN 

END 
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Correspondence with Measured Data 


Ohio State University has, over the past two decades, conducted an extensive 
program of differential cross section measurement for various surfaces. In 
order to verify the theoretical solution, data from this study was compared 
to theoretical predictions. 


The complex permittivity of concrete and asphalt were measured by Peake 
(Ref. 62) by sawing a section of paving material and measuring the 
permittivity in a waveguide bridge. The values he found were: 


Concrete, X-band 
K-band 

Asphalt, X-band 
K-band 


e = 6.5 + j 1.5 
e = 5.5 + j .5 

e = 4.3 + j .1 
e = 2.5 + j .6 


Using these values of permittivity, several comparison runs were made. 
Figures C-6 and C-7 show the best matches found for concrete at X-band. 

Both vertical and horizontal polarization show significant deviations only 
at near vertical incidence. This is not an error in the theory, but is 
merely a result of ignoring the specular component of the reflection. At 
near vertical incidence, concrete at X-band appears quite smooth, and 
thus has a significant specular return. Since Peakes equipment only 
measured over a 1 foot square, the specular component can easily be signifi- 
cant at angles relatively far from vertical incidence. In both horizontal 
and vertical polarization, the best match was obtain for the same roughness, 
.024, but with slightly different slopes. 
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Figure C-6 — Best Match for Backscatter from Concrete at X-band, vertical polarization 
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Figure C-7-- BEST MATCH FOR BACKSCATTER FROM CONCRETE AT X-BAND, HORIZONTAL POLARIZATION 
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Taking the surface parameters which yield the best data at X-band and 
directly scaling them to K-band yield figures C-8 and C-9. These figures 
show a relatively good match between the theoretical prediction and the 
measured data. A few other values for the slope and roughness parameters 
have been attempted, with poorer results. 

Figure C-10 shows the result of a computation for vertical polarised X-band 
backscatter off an asphalt surface. Again, the only significant deviation 
between theory and measurement is at angles which could have significant 
specular contributions. 

APPLICATION OF SCATTERING MODELS TO ILM 

Information on the scattering properties of natural surface is required in 
several areas of the ILM sensor analysis. Scattering data provides a 
means for 

o Differential cross section estimation 
o Multipath environment definition 
o Mitual interference analysis 

Extension of Cross Section Data 

The most obvious usage of scattering models is to extrapolate available 
radar cross section data. Nearly all of the measured cross section data 
available is for incidence angles between 10 and 80 degrees. For the ILM 
program, 84 to 88 degrees are typical incidence angles. Thus at incidence 
angles of interest, very little data is available. Further, the beamwidth 
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Figure C-8 — Backscatter from concrete at K fl band, vertical polarization 
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Figure C-9 — Backscatter from concrete at K band, horizontal polarization 
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Figure C-10— Backscatter franAsphalt at X-band, vertical polarization 


C-36 


aii lay i a a a a a a a a a a i a 



and sidelobe levels of the measuring antenna are not normally known. At 
high incidence angles where the differential cross section is relatively 
high and changing slowly with angle these items are fairly unimportant. 
However, at low grazing angles, the cross section is low and varies rapidly 
with angle. Therefore, antennas of different beamwidth (or systems with 
diffeerent pulse width for a pulse measuring system) can have large 
discrepancies when measuring the same surface, since the major power contri- 
bution may come from slightly different angles. Also, when the cross section 
is very low and the measuring system uses a CW technique, significant 
error contributions can be made by energy entering the sidelobes. Thus, 
even that data which is available at these low grazing angles is of question- 
able validity. Since a computer program is available which provides cross 
section data in good agreement with the high angle measured data, which is 
in some sense mathematically reasonable, and which is intuitively correct, 
the outputs from this program can be used with some confidence. 


The same caveat applies to using the theoretical results as was mentioned 
in criticizing the measurements. The cross section is changing rapidly with 
incidence angle, and thus some gain function must be included in the 
integration of the radar range equation 


P = 
r 



C* fit G R Y dS 

Js r 4 


(C-41) 
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Multipath Analysis 


Traditional multipath analyses have assumed a microscopically smooth, 
perfectly conducting surface. Under this conditions, the energy incident 
on the surface is completely reflected in the specular direction, with a 
phase reversal for horizontal polarization and phase modification for vertical 
polarization. The reflected wave is perfectly coherent, and hence the signal 
at the receiving antenna is purely a function of geometry. 


From the foregoing analysis, it can be seem that although specular multipath 
exists, it is only one component of the multipath. Diffuse multipath 
or that energy received from reflections at angles other than the specular 
is another important contributor. The data on scattering at arbitrary 
angles allows calculation of this effect. 


Since scattering from natural surfaces covers the entire hemisphere, it is 
feasible that some contribution to the diffuse multipath comes from the 
entire area contained by the horizon. In order to simplify the analysis, 
it is necessary to limit the area of multipath contribution. If the very 
rough gaussianly correlated model is used, the mean square value of the 


reflection coefficient for an differential area dS can be represented by: 

cot g o exnf - tan g xv ) (C-42) 

tan 2 $ 0 


<P 2 > 


dS 


4 irr 


4 F 

COS e. 


xy 


where r is the distance from the area to the receiver 
-1 2 <* 

B q =tan — is the mean square value of the slope 

of the irregularities 

and is the angle between the bisector of the incident and 

reflected rays and the Z-axis. 
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In this model, all significant contributions to the received energy come 
from areas where £^<80 . By extending this approximation to all 

models the glistening surface, i.e., that area having significant contri- 
bution to the diffuse multipath signal, may be defined. 

It can be derived geometrically that if the area is not limited by the 
antenna patterns of the transmitting and receiving antennas, the glistening 
surface is approximately a trapazoid with sides defined by (ref. C-l) . 

8=6 

o 

and ends defined by 
d 1 =h 1 cot 2 g o i=r,t 

where d^ is the distance to the start of the glistening 
surface from the ith antenna 

as shown in Figure C-ll. 


Limiting the area of integration to this surface, the mean value of 
received diffuse multipath power is given by the bistatic radar range 
equation: 


P 

r 


V 


(4tt) 


Js Y (6l ’ 92> 


0) VW 2 > °K <vvf > ds 


R 1 2R 2 2 


where 0 t is the transmitting depression angle 


0 r is the receiving elevation angle 


(C-43) 
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Many analyses (ref C-l chapter 7 for example) have shown that the distribution 
of this power is approximately uniform in phase and Rayleigh in power. 

Since any other assumption yields an extremely complex expression depending 
on surface parameters, polarization, etc., the Rayleigh distribution will 
be assumed here. Thus the diffuse multipath appears at the receiver as an 
additive noise term, which provides an upper limit on the signal to 
noise ratio. 

Since the width of the glistening surface is very narrow, the azimuth antenna 
gain variation and change in cross section can be approximated as a 
constant. The width of the glistening surface is: 

w = [ V + £ h 2 tan0 2 (h r h 2 ) ] tan 6 o 

where h^ is the height of the transmitting antenna 

h£ is the height <£ the receiving antenna 

£ is the ground range from the transmitter to the 

receiver 


The total diffuse multipath power is then: 


p i 2 >» h i cotB h 

_ P A r l o —i o—7- -i r -i 

P r = * / Y (tan A ), tan A ( ^ ) , 0) G (tan A ( 0 > 

fA-n\ 3 J 1 2 

k } h_cotB 2 2 

2 ° (t -r) r 


-1 h 2 

G (tan A ( — ) -0 , 0) (rh, + (£-r)h_) 2 tan 6 

r r r 1 2 o 


dr 
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This expression is directly usable for a CW system or a pulse system with 
wide pulses, however, since the path length is varying throughout the range 
of the dummy variable r , this power is spread over time. For a ground 
radar at nearly zero antenna height, the maximum range difference occurs 
at the end of the glistening surface nearest the aircraft. Assuming 
the glistening surface to extend all the way to the aircraft, the maximum 
range difference is from .02 to .1 times the slant range. Thus the maximum 
time difference is from .06 1 to .3 l g microseconds, with lg in Km. In 
a high resolution pulse radar, the integral must be broken up into differ- 
ential ranges, and the diffuse multipath from a given resolution cell 
applied to the appropriate cell. The total multipath is then the sum 
of all the multipath elements at that time difference. 

Mutual Interference Analysis 

One source of mutual interference between imaging radars is the probability 
that the backscattering of a pulse from one aircraft illuminates the 
antenna of another aircraft, giving a phantom image. 

If only the mainlobes of the antennas are considered, this problem is 
very similar to the simple radar problem. However, the general scattering 
equations show that high angle diffuse scattering or specular scattering 
entering the back or side lobes may be of a higher level than direct 
backscatter. Thus, the general scattering equations should be used in any 
mutual interference analysis. 
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APPENDIX D 

ATMOSPHERIC ATTENUATION OF MICROWAVES 


GASEOUS ABSORPTION 

The absorption of microwave energy by atmospheric gases is due to cyclotron 
resonance of the molecules of the constituent gases. In the frequency 
range of microwave sensors, the absorption is due to the 1.35 cm resonance 
of water vapor and a series of resonances of oxygen centered about .5 cm. 


The general theory of gaseous absorption of microwaves has been formulated 
by VanVleck (Reference D-l) and the constants in the formula measured by 
Birnbaum and Maryott (Reference D-2) , Artman and Gordon (Reference D-3) , and 
Becker and Autler (Reference D-4). 


The absorption by oxygen is given by: 


a 

o 



<f4 2 


AV 


1 


2 2 
J/A + AVj 


+ 


AV, 


(2+1 /A ) +AV 


2 

2 


+ 


&V 2 1 
(2-l/A) 2 +AV 2 J 


(D-l) 


where 


a 


o 


A 

P 

T 

AV 

AV 


1 

2 


is the attenuation in dB/Km 
is the wavelength in cm 


is the pressure in atmospheres 


is the temperature in degrees Kelvin 
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The absorption by water vapor is 

5/2 


0) 


. 0318P ( 293 ) 


-644/T 


X 2 


r av 3 

L( 1A --i- ) 

1.35 ; 


2 2 
+ av 3 


(D-2) 


.1/A -t- 1 . 2 2 

( T3 5 ) + 


+ *° 5 V 3 


C^> 


where 


a 

(0 



is the water vapor absorption in dB/Kin 
318 !/ 2 

.087 P (— ) (l+.0046p ) 


P is the absolute humidity in gms/m' 


Many experiments have shown these equations to accurately represent the 
gaseous absorption of the atmosphere. By integrating these functions over 
the transmission path, the total path gaseous absorption may be found. 

A simple model for gaseous absorption in the atmosphere which corresponds 
fairly accurately to observations is the bi-exponential model. In this 
model, the attenuation coefficient due to gaseous absorption at 35 GHz is: 


a a , e 
g = d 


h(68. 6 - 2.75k) 

T a) 


- lOOh 
e 2090 - T 

o 


(D-3) 
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where 


a sf is the ground level dry absorption coefficient 

T 0 

h is the altitude (Km) 

k is the temperature/lapse rate 

T 0 is the ground level temperature (°K) 

T -3 

= (24 - o) Pw x 10 is the ground level wet 

20 101.3 absorption coefficient 

P is the atmospheric pressure (KPa) 

W is the absolute humidity (g/m ) 

This model is useful because of its analytical tractability. Since the 
maximum gaseous attenuation of interest is less than 5 dB, any errors 
caused by the use of this formula will be insignificant. 

ATTENUATION BY CLOUDS OR FOG 

Attenuation of microwave by clouds or fog is of a considerably different 
nature than that due to either rain or water vapor. This is due to the 
scattering characteristics of the very small (<.01 cm) diameter drops. 

Fog attenuation was derived by Gunn and East (Reference D-5) with the results 
shown in Table D-l. Attenuation at frequencies below X band are not 
significant over the path lengths considered, for example at 3 GHz with a 
30 m visibility, it requires a 50 Km path to obtain 1 dB of path atten- 
uation. 

The data in Table D-l is presented in terms of dB/Km/g/m which requires 
knowledge of the amount of condensed water. Based on attenuation measure- 
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merits at several frequencies, an empirical relationship between RVR and 
water content accurate to about ±107» is 

3.04 - 1.37 logV = logW (D-4) 


where 

V is the RVR in feet 
W is the water content in gm/m 

In this relationship, zero visibility would require infinite water content. 

3 

However, the highest reported water content of clouds is about 4.0 gm/m 
in isolated cases of cumulus congestus clouds. Thus, zero visibility fog will 
be considered to have a water content of 4.0 gm/m . Under these assumptions, 
the fog water content of Table D-2 will be used for the analysis. 

Table D-l. Specific Attenuation dB/Km per g/m^ 

Condensed Water in Cloud or Fog 


Temperature, C 

C Band 

X Band 

Ku Band 

Ka Band 


6 GHz 

9 GHz 

15 GHz 

35 GHz 

20° 

.018 

.055 

.111 

.523 

10° 

.024 

.073 

.155 

.586 

0° 

.035 

.101 

.228 

.827 

-8° 

.045 

.131 

.291 

1.043 


D-4 

II I I I I I I I I I I I M M'S I I 1 1 



Table D-2. Assumed Hater Content of Fog 


RVR, ft 


1,200 

700 

150 

0 


g/in of Water 
.065 
.14 
1.1 
4.0 


ATTENUATION BY RAIN 

The attenuation of microwaves by rain is the most significant and 
simultaneously the least predictable of all atmospheric degradations. 

The theoretical foundation for predicting rain attenuation is the paper 
of Ryde and Ryde (Reference D-6) , which assumed a particular distribution 
of water drop sizes and derived expression for attenuation based on Mie 
scattering. The resultant attenuation values can be quite closely 
approximated by a function of the form 

A = kaR b (D-5) 

where 

A is the specific attenuation (dB/Km) 
a is a function of frequency 

k is a function of temperature and frequency 
R is the rainfall rate (mm/hr) 
b is a function of frequency 
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Values of the parameters k,a, and b are given in Table D-3. 


Table D-3. Factors in the Specific Attenuation 
Equation for Rain 


Band 

Coefficient a 

Exponent b 

Temperature Factor k 
0° 10° 18° 

30 

C 6 GHz 

.0027 . 

1.25 

.85 

.93 

1.0 

• 

00 

O' 

X 9 GHz 

.0079 

1.28 

.76 

.90 

1.0 

• 

VO 

O 

Ku 15 GHz 

.033 

1.16 

• 

00 
t— 1 

.92 

1.0 

.93 

Ka 35 GHz 

.17 

1.05 

.86 

.96 

1.0 

.95 


Several measurement programs, summarized by Medhurst (Reference D-7) seemed 
to show a wide variance from this theory. Medhurst rederived the theoretical 
attenuations based on drop size distributions which would yield minimum and 
maximum attenuation limits, and shows that measured attenuations lie 
outside even these bounds. However, the assumption of uniform rain over 
the propagation path is maintained throughout Medhurst's work. In certain 
types of rain, particularly cumuliform rain, there exist large gradients of 
rain fall rate in both space and time. Recent experiments (Reference D-8 and 
D-9) designed to account for these variations have shown extremely good 
agreement with Ryde's theory. In Reference D-8, attenuation of signals 
transmitted from an aircraft to a ground station was correlated with radar 
reflectivity of the rain, using the approximation 


where 


Z = 200 R 1 * 6 


(D-6) 


Z is the radar reflectivity factor. 
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The only attenuation measurements which differed from Rydes theory via 
the reflectivity approximation were measurement where there was evidence 
of hail or snow mixed in the rain. In these cases, the attenuation was 
significantly less than it would be for pure rain, as theory predicts. 

Therefore, it appears that inconsistancies in measured data are due more 
to inaccuracies in measuring the spatial and temporal variations of the 
rain than to any basic fault in Rydes theory. Medhurst's minimum and 
maximum limits are not reasonable for radar performance calculation, 
since it is highly improbably that any rain shower would consist of 
uniformly sized drops, pathologically sized to provide the highest or 
lowest possible attenuation. Any arbitrary variance in attenuation would 
be as mathematically viable as Medhurst's minimum and maximum. 

ATTENUATION BY ICE AND SNOW 

Solid water can be present in the atmosphere in several forms: ice fog 

or cloud, hail, or snow. Because of the different dielectric constants 
of solid water, its attenuation is generally insignificant. However, 
if the solid water is coated with a layer of liquid water, it's attenuation 
can be as great or even greater than the attenuation of the equivalent 
liquid water particle. 

Therefore, the specific attenuation of ice fog or ice clouds will be 
considered to be zero, as will the attenuation of hail or snow if the 
air temperature surrounding the hail or snow is below 0°C. If the 
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surrounding air is above 0°C, the specific attenuation will be considered to 
be anywhere between zero and the value for the equivalent amount of rain. 

ILM WEATHER DATA 

There are four basic climatological conditions to be investigated for ILM. 

They are : 

o CASE 1: Summer Rain, Maritime Tropical Climate 

o CASE 2: Radiation or Advection Fog, Temperate Climate 

o CASE 3: Inland Evaporation - Fog, Temperate Climate 

o CASE 4: Winter Snow, Temperate Climate 

Variations in RVR for cases 2 and 3 will be accounted for by using subcases, 
where 

o Subcase 2 (or 3).l is 1200 ft RVR 

o Subcase 2.2 is 700 ft RVR 

o Subcase 2.3 is 150 ft RVR 

o Subcase 2.4 is 0 ft RVR 

Case 1 is representative of a summer thunderstorm on the gulf coast, with cloud 
tops to 50,000 ft and 16 mm/hr rainfall. Case 2 is representative of coastal 
fog, or fog associated with high pressure cells in the Midwest. Case 3 
is frontal fog, usually experienced in the East and South during the spring. 
Case 4 is a typical winter snow storm caused by maritime polar air over- 
running modified continental polar air. 
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The vertical profiles for the various cases is shown in Figure D-l. 

Based on the weather cases defined above, specific attenuation profiles for 
C, X, Ku and Ka bands have been computed as shown in Tables D-4 through 
D-10. Gaseous absorption was computed using the VanVleck equations. Rain, 
cloud, and fog attenuations were computed by interpolating values from 
Tables D-l through D-3. 

Since the fog for weather case 2 is only 60m thick, it can be assumed to be 
at constant temperature. Therefore, Table D-5 is only the gaseous attenuation. 
The specific attenuation of fog must be added in the first 60m to compute the 
total specific attenuation. Values of fog attenuation for weather case 2 
are given in Table D-ll on page D-15. 
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Temperature in degrees Kelvin 
Hum is absolute humidity in gm/m 











Table D-4 — SPECIFIC ATTENUATION DB^KM* ILM WEATHER CASE 1 


ALTITUDE C BAND X-BANB KU BAND KA BAND 


ILOMETERt 

6 GHZ 

9 GHZ 


1 5 GHZ 

35 GH 

7 

TEMP 'K 

0. 00 

0. 9000E-01 

0.2775E 

00 

0.8412E 

00 

0.3196E 

01 

295. 0000 

0. 10 

0. 9043E-01 

0.2785E 

00 

0.8430E 

00 

0.3199E 

01 

294.3500 

0.20 

0. 9037E-01 

0.2796E 

00 

0.3449E 

00 

0. 3201 E 

01 

293.7000 

0.30 

0. 91 3 IE- 01 

0. 2807E 

00 

0.8469E 

00 

0. 3204E 

01 

293. 0500 

0.40 

0. 9176E-01 

0.2317E 

00 

0.8488E 

00 

0. 3207E 

01 

292.4000 

0. 50 

0.1016E 00 

0.3115E 

00 

0. 9 09 IE 

00 

0. 3475E 

01 

291 . 7500 

0.60 

G.1022E 00 

0.3132E 

00 

0. 9126E 

00 

0.348 IE 

01 

291 . 0908 

0.70 

0.1019E 00 

0.3117E 

00 

0. 9086E 

on 

0. 3468E 

01 

290. 4316 

0. SO 

0.1015E 00 

0.3097E 

00 

0.9035E 

00 

0.3453E 

01 

. 289.7725 

0.90 

0.1010E 00 

0.3078E 

00 

0.3984E 

00 

0. 3438E 

01 

289. 1133 

1 . 00 

0.1005E 00 

0. 3053E 

00 

0.3933E 

0 0 

0.3424E 

01 

288.4541 

1.10 

0.1001E 00 

0. 3039E 

00 

0.3882E 

00 

0. 3409E 

01 

287.7949 

1 .£0 

0. 9963E-01 

0. 3019E 

00 

0. 3832E 

0 0 

0. 3394E 

01 

287. 1357 

1.30 

0. 9918E-01 

0.30 0 OE 

00 

0.8782E 

0 0 

0. 3 33 OE 

01 

286 . 4766 

1.4 0 

0. 9S74E-01 

0. 298 IE 

00 

0.8732E 

0 0 

0. 3365E 

01 

285.8174 

1.50 

0. 98302-01 

0. 2962E 

0 0 

0 . 36S2E 

00 

0. 335 IE 

01 

285. 1582 

1.60 

0. 9739E-01 

0. 2944E 

00 

0.3638E 

0 0 

0.3338E 

01 

284.5512 

1.70 

0.9750E-01 

0.2928E 

00 

0.8596E 

0 0 

0. 3326E 

01 

283 . 96 06 

1.30 

0. 971 IE- 01 

0.291 IE 

0 0 

0.S554E 

00 

0.3314E 

01 

283.3701 

1 . 9 0 

0. 9680E-01 

0.2894E 

00 

0.3513E 

00 

0. 3302E 

01 

282.7795 

9. 00 

0.9659E-01 

0.2876E 

00 

0.8475E 

00 

0. 329 OE 

01 

282. 1890 

3.00 

0. 9470E-01 

0.2699E 

00 

0.8100E 

00 

0. 3172E 

01 

276.2835 

4 . 0 0 

0 . 9492E-0 1 

0. 2686E 

00 

0 . 7938E 

00 

0. 31 19E 

01 

270. 9718 

5. 0 0 

0. 9327E-01 

0. 2693E 

00 

0.8045E 

0 0 

0. 3147E 

01 

265.5238 

6. 0 0 

0. 2033E-01 

0.5701E- 

01 

0. 1258E 

0 0 

0. 4496E 

00 

259 .5714 

7. 00 

0. 1654E-01 

0 . 46' 3 3E - 

01 

0 . 1 02 1 E 

00 

0.365 IE 

00 

253. 0682 

S. 0 0 

0.3352E-02 ' 

0 . 8 3 36E- 

02 

0. 1778E- 

01 

0.63UE- 

01 

246.2572 

9. 00 

0. 2663E-02 

0. 6629E- 

02 

0,141 OE- 

01 

0. 5007E- 

01 

239. 0394 
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Table D-5 


-- SPECIFIC ATTENUATION HB'KMj ILM WERT HER CASE 2 


TITUDE 

C BAND 

X-BAND 

KU BAND 

KA BAND 


LDMETER 

S €• liHZ 

9 3HZ 

15 GHZ 

35 GHZ 

TEMP 'Y, 

0. 0 0 

0. 7765E-G2 

0. 1 04 IE- 01 

0. 2143E-01 

0. 6914E— 01 

282. 7000 

0 . 1 0 

0. 7563E-02 

0. 1 01 IE- 01 

0. 207 0E- 01 

0. 6675E- 01 

282. 0000 

0 . 2 0 

0. 7366E-02 

0 . 982 1 E — 0 

0. 1999E-01 

0 . 644 OE- 0 1 

281 . 3000 

0 . 3 0 

0. 7174E-02 

0. 9533E-02 

0. 1929E-01 

|J • 621 OE— 0 1 

280.6000 

0 . 4 0 

0 . 6937E-02 

0. 9261E-02 

0. 136 IE- 01 

0. 5935E-01 

279. 9000 

0 . 5 0 

0. 6 3 04 E- 03 

0. 399 IE- 02 

0. 1 794E-01 

0. 5765E-01 

279. 2000 

0 . 6 0 

0. 6625E-02 

0. 3729E— 02 

0. 1729E-01 

0.555 IE- 01 

273. 5164 

0 . 7 0 

0.645 IE- 02 

0. 3472E-02 

0. 1665E-01 

0.534 IE- 01 

277.8328 

0 . 3 0 

0. 623 IE — 02 

0 . 8 8 8 E — 0 

m .16' u3E — u 1 

0. 5137E-01 

277. 1492 

6.9 0 

0 . 6 1 1 6E - 02 

0. 7978E-02 

0. 1542E-01 

0 . 4936E-0 1 

276. 4656 

1 . 0 0 

0. 5955E-02 

0. 774 IE- 02 

0. 1433E-01 

0. 474 OE- 01 

£75.7320 

1.10 

0. 5798E-02 

0. 7509E-02 

0. 1424E-01 

0. 4549E-01 

£75. 0934 

1.20 

0. 5644E-02 

0. 7233E-02 

0. 1367E-01 

0. 436 IE- 01 

£74.4148 

1 . 3 0 

0.5495E-02 

0. 7063E— 02 

0. 131 IE-01 

0. 4173E-01 

£73.7313 

1.40 

0. 5349E— 02 

0.6348E-02 

0. 1257E-01 

0. 4 000E-01 

273. 0477 

1.50 

0. 5203E-02 

0 . 6639E— 02 

0. 1204E-01 

0 ■ 3 3 2 5 E — 0 1 

272.3641 

1 . 6 0 

0. 5072E-02 

0. 6449E— 02 

0. 1 16 0E-01 

0. 3682E-01 

271 . 6913 

1.70 

0. 4942E-02 

0. 6269E-02 

0. 112 OE— 0 1 

0. 3552E-01 

271 . 0220 

1 . 3 0 

0. 4314E-02 

0 • 6 093E- 02 

0. 1 030E-01 

0. 3424E-01 

870. 3523 

1 . 9 0 

0. 469 0E-02 

0. 5922E-02 

0. 1 042E-01 

0. 3299E-01 

269. 6835 

3 . 00 

0. 4569E-02 

0. 5755E-02 

0. 1 004E-01 

0. 3177E-01 

269. 0142 

3. 00 

0. 3514E-02 

0. 4293E-02 

0. 6729E— 02 

0.21 06E-01 

262.3213 

4 . 0 0 

0. 2737E-02 

0. 333 0E- 02 

0. 5026E-02 

0. 157 IE- 01 

255.5751 

5. 0 0 

0. 2140E-02 

0. 2592E— 02 

0. 3769E-02 

0. 1 178E-01 

248.7470 

6. 00 

0 • 1 6*3 SE — 0 c! 

0. 2039E-02 

0. 2892E-02 

0 . 9 06 OE- 02 

241. 6736' 

7. 00 

0. 1 323E-02 

0. 1616E-02 

0. 227 OE- 02 

0. 7149E-02 

234.7717 

3. 00 

0. 105 IE- 02 

0. 1283E-02 

0. 1303E-02 

0. 5715E-02 

227.3819 

9. 00 

0. 3357E-03 

0. 1 034E-02 

0. 1457E-02 

0. 4654E-02 

220.9921 
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Table D-6 - 

- SPECIFIC 

ATTENUATION DB 

vKM.ILM WEATHER CASE 

3 . 1 


ALTITUDE 

C BAND 

X-BAND 

KU BAND 

KA BAND 


KILOMETER 

S 6 GHZ 

9 bHZ 

15 GHZ 

35 GHZ 

TEMF "K 

0. 00 

0. 9346E-02 

0. 152 IE— 01 

0. 3165E-01 

0. 1 077E 

0 0 

232. 7000 

0. 10 

0. 9135E-02 

0. 1503E-01 

0. 3126E-01 

0 . 1 064E 

00 

232. 0300 

0.20 

0. 903 OE- 02 

0. 1485E-01 

0. 3038E-01 

0. 1 052E 

0 0 

28 1.46 ij 0 

0 . 3 0 

0. 8379E-02 

0. 1463E-01 

0. 3 05 IE- 01 

0. 1 04 OE 

0 0 

23 0 . 84 0 0 

0.40 

0. 3733E-02 

0. 1452E-01 

0. 3016E-01 

0. 1 02SE 

00 

23 0 . 22 0 0 

0.50 

0. 1234E-01 

0. 255 IE- 01 

0. 541 OE- 01 

0. 1919E 

0 0 

279 . 6 0 0 0 

0 . 6 0 

0. 1S33E-01 

0. 2567E-01 

0. 5457E-01 

0. 1935E 

0 0 

278.3573 

0. 70 

0. 1S32E-01 

0. 2534E-01 

0. 5506E-01 

0. 1 95 1 E 

0 0 

273. 1156 

0.30 

0. 1S32E-01 

0. 2602E-01 

0. 5556E-01 

0. 196SE 

0 0 

2 7 »■" . 7 -• 4 

0.9 0 

0. 1232E-01 

0. 2620E-01 

0. 5607E— 01 

0. 1935E 

00 

276.6312 

1. 00 

0. 1S33E-01 

0. 2639E-01 

0. 5659E-01 

0.2002E 

00 

275.3391 

1.10 

0. 1 $34 E- 0 1 

0. 2658E-Q1 

0. 5713E-01 

0. 2020 E 

00 

2 75 . 1 469 

1.20 

0. 1S36E-01 

0. 267SE-01 

0 . 5763E— 0 1 

0.2038E 

0 0 

3 r' 4 . 4 i.i 4 7 

1.30 

0. 1S37E-01 

0. 2693E-01 

0. 5323E-01 

0. 2057E 

0 0 

273 . 6625 

1.40 

0. 1S40E-01 

0. 272 IE- 01 

0. 538 IE- 01 

0. 2076E 

0 0 

2 r 2 . 92 0 3 

1.50 

0. 1$45E-01 

0. 2756E-01 

0. 5948E-01 

0 . 2 1 0 OE 

0 0 

2 ( 2.17 3 1 

1.60 

0. 1S41E-01 

0. 2769E-01 

0. 5975E-01 

0.211 OE 

00 

2 r" 1 . i 5 0 <' 

1.70 

0. 1S35E-01 

0. 2775E-01 

0.599 IE- 01 

0.211 6E 

0 0 

£71.4226 

1.30 

0. 1S29E-01 

0.278 IE- 01 

0. 6007E-01 

0.2122E 

0 0 

271 . 0945 

1.90 

0. 1223E-01 

0. 2738E— 01 

0. 6024E-01 

0.2129E 

0 0 

2 7 U . 7 6 6> 4 

2. 00 

0. 1S13E-01 

0. 2796E-01 

0. 6042E-01 

0.2136E 

00 

270 . 4333 

3 . 0 0 

0. U37E-01 

0. 2393E-01 

0. 6263E-01 

0.2219E 

00 

267. 1575 

4. 00 

0. U33E-01 

0.233 IE- 01 

0. 6253E-01 

0.2213E 

00 

264 . 2 066 

5. 00 

0. 1S69E-02 

0.2435E-02 

0. 4134E-02 

0. 1332E- 

01 

253. 1250 

6 . 0 0 

0. 167 IE- 02 

0. 2 03 IE- 02 

0. 2924E-02 

0 .91 79E- 

02 

242. 5000 

7.00 

0. 1 S23E-02 

0. 1616E-02 

0 . 227 OE— 02 

0.7149E- 

02 

234. 7717 

3. 0 0 

0. 105 IE- 02 

0. 1238E-02 

0. 18036-02 

0. 5715E- 

02 

227.3319 

9 . 0 0 

0. 8357E-03 

0. 1 034E-02 

0. 1 457E-02 

U . 4 6 5 4 E — 

02 

220.9921 
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Table D-7 


— SPECIFIC ATTENUATION DE-KM* ILM WEATHER CASE 3.2 


LTITUDE 

C BAND 

X-BAND 

KU BAND 

KA BAND 


I LOME TER 

S 6 GHZ 

9 bHZ 

15 GHZ 

35 GH 

Z 

TEMP 'Y. 

0. 00 

0. 1 1 17E-01 

0. 2075E-01 

0. 4344E- 01 

0 . 1 522E 

0 0 

282.7000 

0 . 1 0 

0.11 06E-01 

0. 2069E-01 

0. 4338E- 01 

0. 1 52 OE 

0 0 

282. 0800 

0 . 2 0 

0. 1 096E-01 

0. 2065E-01 

0. 4334E-01 

0. 1519E 

00 

281. 4600 

0 . 3 0 

0. 1 086E-01 

0.2 06 IE- 01 

0. 4332E-01 

0. 1518E 

0 0 

280. 8400 

0 . 4 0 

0 . 1 076E- 0 1 

0 . 2 058E- 0 1 

0.433 IE- 01 

0. 1518E 

0 0 

280. 2200 

0 . 5 0 

0. 1234E— 01 ' 

0. 255 IE- 01 

0. 541 0E-01 

0. 1919E 

0 0 

279.6000 

0 .60 

0. 1233E-01 

0. 2567E-01 

0. 5457E-01 

0. 1935E 

00 

278.8578 

0. 7 0 

0. 1232E-01 

0. 2584E-01 

0. 5506E-01 

0. 195 IE 

00 

278. 1156 


Table D-8 

-- SPECIFIC 

ATTENUATION 

DE'KMjILM weather case 

3.3 


ALTITUDE 

C BAND 

X-BAND 

KU BAND 

KA BAND 


KILOMETER 

S 6 GHZ 

9 GHZ 

15 GHZ 

35 GH 

y 

TEMP 'K 

0. 0 0 

0. 3453E-01 

0. 9163E-01 

0. 1943E 00 

0. 7217E 

0 0 

282.7000 

0 . 1 0 

0. 35 07E-01 

0 . 9325E- 0 1 

0. 1986E 00 

0.7359E 

00 

28 2. 0800 

0 . 2 0 

0. 3562E-01 

0. 9487E-01 

0.2029E 00 

0.750 IE 

00 

28 1 . 46 0 0 

0 . 3 0 

0 . 36 1 8E- 0 1 

0. 965 OE- 01 

0.2073E 00 

0.7644E 

00 

230. 8400 

0 . 4 0 

0. 3674E-01 

0.98 13c- 01 

0. 21 16E 00 

0. 7787E 

00 

280. 2200 

0 . 5 0 

0. 1234E-01 

0.255 IE- 01 

0. 541 0E-01 

0. 1919E 

00 

279. 6000 

0 . 6 0 

0. 1233E-01 

0. 2567E-01 

0. 5457E-01 

0. 1935E 

00 

278. 8578 

0 . 7 0 

0. 1 232'E-0 1 

0 . 2534E-0 1 

0. 5506E-01 

0. 1951E 

0 0 

278 . 1 1 56 


Table D-9 --SPECIFIC ATTENUATION DB-'KM. ILM WEATHER CASE 3.4 


ALTITUDE C BAND X-BANP 

KILOMETERS 6 GHZ 3 GHZ 

0. 00 0.1 05 IE 00 0.3053E 00 

0.10 0.1076E 00 0.3124E 00 

0 . 5 0 0 . 1 1 0 1 E 0 0 0 . 3 19 1 E 0 0 

0.30 0.1127E 00 0.3257E 00 

0.40 0.1152E 00 0.3324E 00 

0 . 5 0 0 . 1 2 34E - 0 1 0. 255 IE- 01 

0 . A 0 0 . 1 233E-0 1 0 . 2567E-0 1 

0 . 7 0 0 . 1 232E- 0 1 0 . 2584E-0 1 


KU BAND Kh BAND 

15 GHZ 35 GHZ TEMP 'K 

0 . 65 02E 0 0 0 . 2442E 0 1 282 .700 0 

0.6676E 00 0.2500E 01 282.0800 

0 . 685 0E 0 0 0 . 2557E 0 1 £E; 1 . 46 0 0 

0 . 7 085E 0 0 0 . 26 1 5E 0 1 28 0 . 84 0 0 

0.7199E 00 0 . 26 ZEE 01 280.2200 

0. 54 1 0E-01 0. 1919E 00 279.6000 

0 . 5457E-0 1 0.1 935E 0 0 278 . 8578 

0. 5506E-01 0. 195 IE 00 278.1156 
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Table D-10 

-- SPECIFIC 

ATTENUATION DB 

-K M- ILM WEATHER CASE 4 


ALTITUDE 

C BAND 

X-BRND 

KU BAND 

KR BAND 


KILDMETEP 

S 6 GHZ 

9 GHZ 

15 GHZ 

35 GHZ 

TEMP '¥■ 

0. 00 

0. 8527E-02 

0. 1 019E-01 

0. 1506E-01 

0. 463 0E- 01 

268. 0000 

0 . 1 0 

0 • yjl4E~ 0 cl 

0 . 993SE-U2 

0. 1465E-01 

0. 4504E-01 

367.3400 

0 . 2 0 

0 . 3 1 07E - 02 

0.9677E-02 . 

0. 1 425E-0 1 

0. 438 OE— 01 

266.6800 

0 . 3 0 

0. 7905E-02 

0.943 IE- 03 

0. 1 336E-01 

0. 426 0E- 01 

366 . 02 0 0 

0 . 4 0 

0. 77 09 F- 0£ 

0.919SE-03 

0. 1343E-01 

0. 4143E-01 

365 . ofc. 0 0 

0 . 5 0 

0.7517E-02 

0. 8953E-02 

0. 131 IE- 01 

0. 4027E-01 

364.7000 

0 . 6 0 

0. 7.308E-02 

0. 8705E-02 

0. 137 IE- 01 

0. 3905E-01 

264. 3387 

0 . 7 0 

0 . 7 1 04E - 0£ 

0. 8459E-0S 

0. 1233E-01 

0.3736E-01 

263. 9773 

0 . 3 0 

0. S906E-02 

0.3230E-02 

0. 1 195E-01 

0 . -■ 6 7 OE — 0 1 

263 . & 1 6 0 

0 . 9 0 

0 .671 4 E — 0 S 

0. 7983E-02 

0. 1 153E-01 

0. 3556E-01 

263. 2547 

1 . 0 0 

0 . 6527E- 03 

0 . 77 6 1 E- US 

0. 1 123E-01 

0. 3446E-01 

262.3934 

1.10 

0. 6345E-02 

0. 754 IE- 03 

0. 1 038E-01 

0. 3339E-01 

262 . 532 0 

1.30 

0 .6-1 68E-02 

0. 7327E-0S 

0. 1 054E-01 

0.3234E-01 

262 .17 07 

1 . 3 0 

0. 5996E-02 

0.71 19E-02 

0. 1 021E-01 

0. 3132E-01 

261 . 8094 

1.40 

0. 5329E-G3 

0. 6917E-03 

0. 9885E-02 

0. 3032E- 01 

261 . 4480 

1 . 5 0 

ij . 5666E- 

0. 6720E-02 

0.957 IE- 02 

0. 3935E-01 

261 . 0867 

1.60 

0. 5495E-02 

0 ■ 6527E - 02 

0. 9344E-02 

0. 2869E-01 

261 . 0499 

1.70 

0.5334E-02 

0. 6339E-G2 

0. 9 149E-02 

0. 2814E-01 

261. 1155 

1 . SO 

0. 5159E-03 

0. 61 57E-0S 

0. S960E-03 

0. 3761 E- 01 

261. 1811 

1 . 9 0 

0. 4999E-03 

0. 5 9 32 E- 02 

0. 3777E-02 

0. 2709E-01 

261.2467 

2. 0 0 

0. 4344E-03 

0.581 IE- 02 

0. 8599E-02 

0. 2659E-01 

261.3123 

3. 0 0 

0. 3542E-02 

0. 4374E-02 

0.7039E-02 

0. 2228E-01 

26 1 . 9635 

4 . 0 0 

0. 3709E-03 

0. 3 3 32 E- 02 

0. 5237E-02 

0. 1646E-01 

256.9718 

5 . 0 0 

0 . 3 095E - 03 

0. 2557E-02 

0. 3325E— 02 

0. 1200E-01 

250. 9702 

6. 0 0 

0. 1 664E-02 

0. 2019E-02 

0 . 2882E- 02 

0. 9039E-02 

242.7357 

7 . 0 0 

0. 1 313E-02 

0. 1604E-02 

0. 3258E-02 

0.71 14E-02 

235.4751 

3. 0 0 

0. 1 044E-02 

0. 1279E-02 

0. 179 IE- 02 

0. 5673E-02 

223.5066 

9 . 0 0 

0. 3262E-03 

0. 1 023E-02 

0. 1442E-02 

0. 4609E-02 

231.9449 


Table D-ll. Specific Attenuation of Fog for Weather Case 2. (Db/kn i) 


Subcase 

RVR 

C Band 

X Band 

Ku Band 

Ka Band 

2.1 

1200 ft 

.0015 

.0047 

.010 

.038 

2.2 

700 ft 

.0033 

.010 

.022 

.082 

2.3 

150 ft 

.026 

.080 

.171 

.645 

2.4 

0 ft 

.096 

.292 

.620 

2.34 
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APPENDIX E 


RADIOMETRY COMPUTER PROGRAMS 

A set of computer programs has been written for use in the analysis 
of microwave radiometry. The main programs are: 
o Skytemp 
o Pathtemp 
o Bnis 

The first program, Skytemp, performs the integration of specific 
attenuations to obtain the radiometric sky temperature at various 
incidence angles. It assumes a flat earth and a layered atmosphere, 
and assumes that all significant attenuation occurs in the first 
lOKm of atmosphere. Thermometric temperatures and specific attenuations 
are read from files pre-stored on the H-6080 disk file system, and sky 
temperatures are output to the disk on a file names FLTSKY in a format 
which is easy for the computer to use in further calculations . 

Pathtemp is a very similar program, which integrates the specific atten- 
uation to obtain the total one way attenuation on any glide path from 
any altitude to the ground. It also computes the path emission observed 
at any altitude (under lOKm) and at depression angles from .01745 rad to 
.157 rad (1° to 9°), for the ILM weather cases. Its output is to a disk 
file named FLTPATH for path temperature or FLATRANS for path attenuation. 
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EMIS uses as inputs a description of a statistical rough surface 
(roughness, slope, and perni.ttivity) , sky temperatures, and ground level 
thermometric temperatures. It computes the integrated sum of same sense 
and cross polarized reflected sky emission, and the diffuse percentage 
of power reflected from the hemisphere for both horizontal and vertical 
polarization based on exponentially correlated random rough surface scattering 
theory. The radiometric temperature contributions of diffuse sky, specular 
sky, and emissivity are computed and summed. The results are output to 
a temporary disk file for later printing. 

Listings of the programs and results of interest are attached. Included 
are: 

o Total oneway path attenuations for 2,6, and 16 Km ranges on 
one to nine degree glide slopes 

o Perceived radiometric path temperature for the same conditions 
o Apparent surface radiometric temperatures for selected glide slopes 
and grass, concrete, and snow surfaces 

Grass is described on the printout as a surface with roughness = .5, correl- 
ation length = 3.2, and permittivity 10+jl0. Concrete has roughness = .084, 
correlation length .2, and permittivity 5.5+j.5. Snow is roughness .1, 
correlation length 1.0, and permittivity 3.2+j85. 
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Program Sky temp 


10 DIMENSION ALPHA (28*10) , TH I K ( 26 ) » THFD < 1 8 ) , TSKY ( 1 8 , 1 0 ) * TEMP (28 , 10) 

030 - DATA THIK/20*. 1 ,e*l./ 

AA CALL ATTACH <20, 'DOOOAfl/FLALPHA ? • , 3,0, I S» ) 

A 5 CALL ATTACH(21 , ‘POOOAS/FLTEMP; • ,3,0# ISTATO 

A 6 REV. 1 MD 20... _ 

A 7 PEWIND21 

A3 READ (20»?A1 ,FMD=A9) ( (ALPHA <IO)«I = 1«28),J=1,10) 

A 9 ...A 9. RFAD(21i2A2,F.ND=50) ( (TE.MF ( I ♦ J) ,1 = 1,28) , J= 1 « 10) . . 

50 50 DO 200 <=1*18 

060 IF (K.C-E. 10) GO TO 100 

.. .070 THEP ( '< ) =K .. . ' . . .. 

0S0 THFTA=K/57.2P 

090 GO TC 120 

.. 100 100 THSPU) =K-9 

110 THETA=THEP(K)/5.729 

120 120 CONTINUE 

• ..130 SINTH-SIN (THETA) :.. 

1 AO DO 190 L=1 .10 

1 A 1 IF(K.NE.1)G0 TO 1A5 

. , ’A5 1A5 CONTINUE . . 

6 ATTFNJ=1. 

■ 1A8 TSKY ( < , L ) =0.0 

150. DO 190 1 = 1 *2.8 . I 

160 ATTFNI=ATTFNJ*EXP(-.23*ALPHA(I,L)*THIK (D/SINTH) 

170 TSKY (< ,L) =TSKY (K ,L) ♦ ( 200 , ♦ TEMP ( I , L ) ) * ( ATTFN J-ATTEN I ) 

180 ATTFNJsATTENI ... .'. 

190 190 CONTINUE 
200 200 CONTINUE 

__ 210 CALL ATTACH ( 22 i .' D000AS/FLT5K Y » , 3 » 0 , I STAT , ) ... 

220 RFWIND22 

2 30 WRITE (22,2A0) ( {TSKYU , J) , 1 = 1 , 18)., J=1 , 10) 

_..2A0 2A0 FORMAT (6E12. A) 

2 A 1 2 A 1 FORMAT (2 ( 10E12.A/) ,8F12.A) 

2A2 2A2 FORMAT (1AFA. 2) 

2 A 6. CALL DETACH (20, IS,) J 

2A7 CALL DET ACH (21,150 

2A8 CALL DETACH ( 22 , I S , ) 

...250 STOP . — ■ - - ... • 

260 END 


ORIGIN AD PAGE IS 
OF POOR QUALITY! 
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Program Pathtemp 


0010 DIMENSION TAUA (28) , PATHT (20) , ALPHA (28) ', TEMP ( 28 > »THIK (28) 

. 0020 .DATA THIK/20*. 1 ,8*1./ ..... 

0030 CALL ATTACH(20, 'D00048/FLALPHA* • ,3,0, IS*) 

OOAO CALL ATTACH(21 , • D0004 8 /FLTEf'P ? • *3*0»IS») 

— .0050 CALL ATTACH ( 22 » •D000A8/FLTPATH; • »3*0i IS*) 

0060 CALL ATTACH(23. «D00048/FLATRANS: • .3.0. ISt) 

0070 REV.' I ND20 

— 0080 REWIND21 .... 

0090 REWIND 22 

0100 REWIND 23 

—.01.10 DO 290 . I W X = 1 » 1 0 . - - 

0120 READ C 20*300) ( ALPHA ( I ) *1 = 1*26)’ 

0130 READ (21 *310) (TEMD(I) ,1=1,28) 

- 01 AO DO 280 IC-S=1 ,9 - 

0150 THETA=IGS/57.29 

0160 CSTH=SIN (THETA) 

-0.1.7 0 ptran=i.o : ;.... 

0180 PTMR=0.0 

. 0190 DO 250 J= 1 * 2 8 

-.0200. . FXPFAC = EXR (-. 2 30 1*ALPHA ( J) *THIK ( J) /CSTH) 

0210 P TP A N = ° TR AN * E X PF AC 

0212 IF (PTRAN.CT.O. ) GO TO 220 

-.0 TAUA (J) =99. 99999 

02 j. 6 GO TO 222 

0220 220 TAUA ( J) ~10.*AL0G1C (PTRAM) 

—0222 .222 CONTINUE • - - 

0230 PTMP= (PTNP-200.-TEMP ( J) ) «EXPFAC*200.*-TEMP ( J) 

02A0 PATHT •< J ) =PTMP 

_.0250 .250. CONTINUE — •• 1- ... - 

0260 WRITE (22,320) (PATHT < I ) *1=1,28) 

0270 WR1TE<23, 330) (TAUA(I) *1=1,28) 

0280 280 CONTINUE -• - : 

0290 290 CONTINUE 

0300 300 FORMAT (2 ( 10E12,A/) , 8E 1 2 . A ) 

D310..310 FORMAT ( 1AFA. 2) 

0320 320 FORMAT ( 5E 1 A , 6 ) 

0330 330 FORMAT (5E1A. 6) 

— 0340 -CALL DETACH ( 20 * I S* ) 

0350 CALL DETACH (21 , I S, ) 

0360 CALL DETACH (22, IS,) 

.—0370 CALL- DETACH (23,1 S,) ■ - 

0380 STOP 

0390 END 
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Program Emis 


—I 5To TTTPensTon GAMSUM<2*18) *REF0UT<18) *GLTMPC10) *EMSV(2) t SPFC<?> «SPSKY<10)_ 

1.1 DIMENSION SPCSK Y (9*10) 

012 DIMENSION TSK Y ( 9 * 10 ) * WXC5 < iO) »TD ( 10* 2 ) 

014 DIMENSION TAPPUO, 2) .TSPECdO. 2) «TEMIT(10*?) .EMIT(2) 

020 DATA RSH«5LP.*PRMPE*PRMIM/»5 *3.2* 10# « 10#/ 

025 DATA GLTMO/295. *8*282.7*268*/ 

26 CALL ATTACH (2 1 * • D00048/PUTPUT ; * * 3 *0* 1 SO 

027 REWIND 21 

030 DATA D1MC/89./ 

31 CALL ATTACHI20* •D00049/FLTSKY? • *3*0* 1ST AT ♦ ) 

32 REWIND 20 

033 READ I 20*34) ( (SPCSK Y ( I ♦ J’> 1 1 =1 .9) . t TS* Y (K . J) ,< = 1 , 9) . J=1 » 1 0) 

034 34 FORMAT (6F12. 4) 

035 C40 PRINT 36 

036 036 FORMAT ( "POUCH SURRAC 17 PARAMETERS •PGH«SLP«PRMR=-,PRmtM" ) 

037 READ 38.* DA T 1 * 0AT2 *'DA T 3 *0AT4 

038 038 FORMAT (4cb. 4) 

039 V 5PPC=90-r INC 

040 a INC = DIN'C/57.29 

041 I F (DA T 1 . ro, 0. ) GO TO 43 I RGH = PA T 1 ? SL ° = DAT2 ? PRMRC = D AT 3 ! PPM 1 M = DA T4 

42 DATA WXCS/l.P*?.l*2.2.2.3»?.4»3.1*3.2»3.3*3.4«4.0/ 

043 043 FNIT2=C. 

044 

045 

050 
060 . 

070 
050 
030 
100 
1 10 
120 
130 
140 
150 
1 60 
170 
180 
190 
2 00 
210 
220 
230 
340 
350 
360 

370 

371 

372c ****V C RY rough fy ddnenT I a,L**** 

373 TFMP1=G**?+39.4748*ELP**?*F.XYSG 

374 d FST = 78.9 c 6*C*S|_p**2/ < (CSIN'C+CSR^F) * *2*TF* PI *SORT < TFMP 1 ) i 

375 CO TO 400 

376C ****SL IC-JTL Y RP'JGM cypo^r^Ti r.L*«** 

377 377 SUMG=0.d 

378 SmlST = 0. 

379 FACT«=l. 

380 GTOM=l. 


RUM= 1 00. 0 

F V I T 1 = 0 . 0 

DO 428 j= 1 , 1° 

GAMSIJM ( 1 , J) =0. 

GAMSUM (2 * J) =0. 

RFFCUT (J) = 10*J-O5 

ACFFL =RFF0UT ( J) /57. ?9 
oo 42? 1=1*3 
AO A 7 = ( I-l) / 5 . 7 29 
CSINC=COS(AP'C) 

CSRFF=CD5 (AR CC L) 

C?OA7=CCS (AOA7) 

SNINC = SIN (A INC) 

SNR FF = S I w ( AR e FL ) 

SN0A.7 = S I N ( A0A7 ) 

A 2= (CSI NC*SMPFF + SNINC*CSRFF*CSOAZ) *«2 
A3= (5NI M C*CS D EE+CSINC*SNR“F*CS0AZ) **2 
A8 = SNINC«8.NREF*CS0A7-CSINC*CSRFF 
C8L IN= ( 1 .-AS) / 2 . 

SML IN= ( 1 .+A8) /2. 

CALI FRA»j r L (SNLIN*CSLIN«PR.MPF,PP M H«PV*RH) 
FSAM=A2*(A3+( SNRE r# SNOAZ) **2) /<(l+A5)**2) 

FOR 5= (SMI M C*SNOAZ) **2* < A3*5NOA7**2 ) / ( ( 1 + A5 > **?) 
FXYSD = SNI mC**2 + SNREF**2-2*SNINC*SNR c 'F*C50AZ 
C= (6. ) «*R0H* (CSI‘'C*CSRFF) } **2 

IF (C.LT.12.5) G° to 377 
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38 1 
3*2 
3*3 
384 
3*5 
386 

397 

388 

389 

390 

396 396 
3976 

398 398 

399 399 

400 400 

401 
410 
420 
430 
440 

450 450 

460 800 

470 

490 

500 

510 

520 420 

53C 

540 

550 425 

560 

570 

5 °0 
590 
670 
672 

674 

675 
680 

6 90 
700 
710 
720 

730 666 

740 

742 

744 

746 

750 668 
752 

754 754 
7556 
756 755 

758 

759 

760 
7706 


00 39P >= 1,20 
GTOM = C-TC"A*G 

FACTM=FACTM*m 

TFMPl=l439.4748*FYYSQ* (SLR/M) **2 
Sl'MG=5U^G*C-T0M/ (FACTM*N1*#2*TFVP1*S0RT (TF.MP1 ) ) 

SU M D I F= A3S ( SUMG-SMLST ) 

IFISUvOIP.LT, .01*SUNG) GO TO 399 
IF(v.F0.20)O9JNT 396 ,R P c OUT ( J 5 , S'JMO T F 
5MLST=SUM0 

FORMAT ("SERIFS CONVERGENCF FPROP , ANGLF OF 9F F= " , F5 . 2 . " SUMO I F = " , 

FI 1 .4) 

CONTINUE 

RP5T = 78.956«Fvp (-G) *SUMG* ( 5LP / ( C S I NO *C SRF r ) ) **2 
contimup 

SLOP=SLP/ ( 1 • 4 1 4*RGH) 

IF (CSINC.GT.CSREF) GO TO 450 
TEMPI =APS (SNPPF/CSRPP) 

Sh ADR = SHAO (SL0P»TFMP1 ) 

GO TC 800 
SSAOR = 1.0 

Tfmpi=a9S (SN^'C/CSINC) 

SHAD1=SHAO(SLOP,TFMP1) 

RrST=PF5T»SMA0r*5 u A0R 

GAMSUv ( 1 , J) r (FSAM*RV+FCRS*RH) *RFST ♦GANSUM I l , J) 

GAM SU V (2,j) s (rsftM#Pt;4FC f> 5*RV) #RFST*G AMSUM (2 , J) 

CONTINUE 

EMIT ( 1 ) =FMIT ( 1 ) +GAMSUM ( 1 , J) 

EMIT (?) =EMIT (2) ♦GAMSUM (?« J) 

CONTINUE 

EMIT ( 1 ) = p m I T ( 1 ) *2.424F-3/CSIMC 
EMIT (2) =EMIT (?) *2.424F-3/CSlNC 

CALL FR 5N r L (SMINC,CSINC*pR.mpE,PPMIM,S^FC ( 1 ) .SPEC (?) ) 

RF ST = FXD (-4*39.47 48* ( PGH*C S I NC ) * * 2 ) 

00 668 L = 1 , 2 

SPEC ( L > =SPEC (L) *PEST 

FM5V(L ) =1.-FMIT (L) -SPEC (L) 

IE(FMCV(L) .LT.C. ) PMSV (L) =0. 

HO 66? K = 1 , 1 0 
T 0(K iL) =0.0 
no 666 J=l,9 

G A MT'OT = G A“ SUM (l«9*J) ♦ GAM.S'.'M (l ,10-J) 

TO(X,L)=Tr(K,L) +GAMTOT*TSK Y ( J,K ) 

CONTI mijf 

TO IK ,L ) =Tn (< «L ) *2.4?4F-?/'"MNC 
TEN IT (K «L) sGLTmd (< ) *FVS'/(L) 

TSPEC (K , L ) =SPFC (L) *SPCSK Y (K^PFC ,K) 

T/»PP (4 ,L ) = TO(k- ,l ) + TEMIT (K .L) ♦ T SDPf (K,L ) 

CONTINUE 

VP 1 TP ( 2 1 ,7 84 ) c?3H • 8LP , pc«i».;,r, I»!C 

FORMAT (/"SURFACE raPA M ET r PS: °0UGHNP8 p = " ,P6. 3 , " COPULATION LE*C-TH = ". 

F5.2," PEP*’ 1 TTVITYr" ,F5. 1 ,?M4j,F4. 1 lMCT r 'FNCF ANGL F= " , F4 . I ) 

FORMAT ( "orovi ti viTY = " ,F4. 1 » ?K* J , F4. 1 , " I \C T OP NC C ANGLF=" ,F4. 1 ) 

WRITE (21 ,825) 

V.'PITF (21,826) 

WP I TF (21 ,695) (WXCS (M ) . (TO (< ,L ) ,TSPPO (V , L ) • TE« I T ( K , I ) , T APP (K ,1 ) , 

L = 1 • 2 ) ,< = 1,10) 
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700 470 FC' d “AT ( "TUFT Af V r o T I r A I_ h'OP T ■ , C'! T ' t " ) 

800 4 8 0 rpo m ft T (F5«l « 4 X , c 1 1 .4»1X.F1 1 • M 

810 695 FpcwftT ( 1 0 <2Y*F3. 1 .a <5X , c 5. 1) ,7X ,4 (SV ,F?. 1) / U 

8?5 825 FORf.ftT ( t h! , c- ATF'E' 5 • ? * ’ FHP- 1 RFU5F . 2 X * 8USPFC UL 4P * 2X . 8MF** t TTEO ,2X, 

8266 FHArraopr T . 1 ox . Ol io TFF'J8 C .?v . 8WSDFCI IL AR . ?X , FURm j TTFT) ,2v»8HAPPARENT) 

827 826 FORMAT (" C A5 E , 6X . 3HSF Y » 7Y , 3HSF Y , 1 6X . "TFNPFR ATUR^" . 1 OX , 3nSK Y , 

8 2 86 7X, 34SXY* 1 5X." Tr. -.ip.cq AT i |dt ••/] OX, "VERTICAL POL AR T 7 AT T ON" , ?0X , 

e29 f . 2 0 X , "HOP 1 7.P.N T AL n 0L AR I 7 AT I Of:" ) 


830 RFAp 4 8? ,0 INC 

840 4P2 FORMAT (F10. 4) 

8 50 IT (OIYC.NF.C. ) F-n TO 40 

851 CALL DFTAFH f 2D « I ST AT ,) 

852 CALL DETACH ( 2 1 , T STA.T , ) 

858 CALL DETACH (2 1 , I S, ) 

860 STOP 


870 END 

880 FUNCTION SHAT ( SLOP , T A VA ) 

990 t NET = 8L0P/TAN'A 

900 E°FC= • 5/ ( ( 1 *.27830 3 *THET + .2"»0389*THFT**2*.000972*THFT**3*. 078108# 

9106 fHFT**4)**4) 

920 SHA0=EXP (-.5*TANA*E op C> 

930 870 PFTIJRN 


940 880 FORMAT ("SHADOWING OUT ST T ONAPLE ) 

950 END 

1000 SUPPOUTINE FPSNEL (SNLIN,CSLIN.,PP.MRE*PR M I>,RV,RH) 

1010 P1=PPMPE-SNLIN 

1020 A6S0=SQRT <P1**2+PRMIM«*2) 

1030 R2=.5#ATAN2 (PRMIM/Bl ) 

1040 CSB=C0S (P2 ) 

1050 SNP = S I N ( p 2 ) 

1060 PNOR Y=PRMRF**2*PRM I M**2 

1070 PV= ( (PN0RM*CSLIN)**2*4*A6SQ* <PRMRE*SNB> **2*CSL IN*A6S0»*?-2* 

10806 A6SG*PN0RM*CSLIN) / ( <PNORM*CSLIN*2* 

1090 6 SORT ( A6SQ*CSL IN)* (PRMRE*CS8*PRMI M#5NB) *A65Q) *#2 ) 

1100 CONTINUE 

1110 RH= <CSLIN**2*46SO*4»SNB**2*C5LIN*A6SQ**2-2*A650*C5! TN> /< <CSl IN* 

11206 2*SQRT< 

1130 6 A6S0*C5L IN) *CSB#A6SQ) *#2 ) 

1140 RETURN 

1150 END 
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APPENDIX F 

MLS CONFIGURATION K AIRBORNE EQUIPMENT 


PHYSICAL DESCRIPTION OF FAA-K AIRBORNE EQUIPMENT 

The H-80 Airborne Equipment Set offers operational and installation flex- 
ibility through compact, modular equipment packaging. 


The H-80 airborne system will meet all FAA-K equipment requirements. Figure F-l 
shows two airborne sets in a typical redundant aircraft installation. 

Features include a standardized package design and simplified interconnections 
that will permit straightforward installation of the equipment in the DC-6 
or CV-880 aircraft. When the airborne equipment is used in the dual con- 
figuration illustrated, provisions are made for interconnection of the 
two HN-700 Angle Receiver/Processors to implement built-in cross monitoring 
capability. 


The physical characteristics of each H-80 subsystem are summarized below and 
discussed in detail in the following paragraphs. 

Physical Characteristics of H-80 Equipment Set 



Angle Rec/Proc 
(HN-700) 

DME Interr. 
(HN-800) 

Control Unit 
(HC-400) 

Antennas 

(HL-181) 


Size 

3/4 ATR Long 

1/2 ATR 
Short 

5.75"W, 
4. 125"H, 
2 ,5"D 

■ 

Not Defined 

It 

»t 

Weight 

20 lbs 

10 lbs 

3 lbs 

1 lb 

Not Defined 

Power 

145 w 

40 w 

15 w 

— 

— 
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Figure F-l - - Redundant Installation 
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Although the HR-800 DME Interrogator is considered an integral subsystem, 
the HN-700 Angle Receiver /Process or can operate independently. 

The H-80 airborne set is designed to use about 200 watts of 400 Hz, single 
phase, 115v aircraft power. The HN-700 connects to the aircraft power source 
and supplies the required dc power to the other units as required. 

The airborne units have been designed for hard-mounting in the aircraft, and 
will operate and withstand expected prototype aircraft environmental conditions: 
Temperatures -40 degrees to +65 degrees C 
Vibration +2 g's per MIL-STD-810B (Curve B, Figure 514-1) 

The HN-700 is packaged in a standard 3/4 ATR long case and has been dimen- 
sioned for hard mounting to the test aircraft equipment rack for ease of 
installation. Mechanical holddown clamps allow quick installationand removal 
of equipment. Plug-in circuit boards allow fast replacement of defective 
circuits for ease of maintenance. 

System interconnections have been simplified and provide ready access to 
connectors and test points for in-flight equipment monitoring during the 
test program. These connectors, as well as the test point access, are 
located on the front panel. The HN-700 interconnects to the HC-500 control 
unit for channel select and azimuth and elevation path selection. Additionally, 
The HN-700 provides basic low dc voltages, RF signal down conversion, and a 
frequency synthesizer signal for the HR-800. The HN-700 angle deviation out- 
puts are scaled to standard ILS course width sensitivities and are therefore 
compatible with the existing set of avionics in the test aircraft. 
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The HR-800 has similar packaging, mounting, and interconnect features and 
is housed in a 1/2 ATR short case. Case size (and cost) has been sub- 
stantially reduced by a sharing of the HN-700 RF components. 

The Control Unit (HC-400 provides a system on-off switch, a system go/no- 
go status light, azimuth and elevation tracking status lights, channel select 
and MLS path select. The HC-400 unit is of modern keyboard type entry with 
gas discharge display characters formed against a black background for 
optimum readability. Polarized non-ref lective filters combined with proper 
contrast ratio and brightness make the display readable in any ambient 
light conditions including direct sunlight. The control unit is sized to 
fit standard instrument panel mounting slots and is held in place with 
Dzus fasteners. 

Three antennas are used for both equipments. In a dual installation both 
H-80 sets may share the same three antennas. 

The HL-181 directional antenna is a conventional horn intended for mounting 
externally on the nose, or inside the aircraft radome. It is about 4x5 
inches with a depth of 6 inches, including the fiberglas cover. 

The HL-362 and HL-363 omnidirectional antennas are not defined and must be 
customized to the aircraft considering available locations, aircraft 
geometry and MLS equipment location. Typically it is expected that "thimble- 
sized" linear stub type antennas will find wide application. 
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PERFORMANCE CHARACTERISTICS OF FAA-K AIRBORNE EQUIPMENT 

The performance characteristics of the configuration K-FAA Airborne Doppler 
MLS Equipment Set, designated H-80, will meet the FAA requirements for high- 
capability guidance equipment for aircraft engaged in autoland operations 
at primary hub airports. 


The H-80 airborne equipment has been designed to provide precise takeoff and 
landing terminal area guidance information, under Category III weather 
conditions, for fixed-wing civil aircraft operating with autoland avionics at 
suitably equipped major runways. The airborne set will provide the accuracies 
and functional characteristics in Tables F-l and F-2. 

The HN-700 provides the 200 channel frequency synthesizer, down conversion 
of both the angle guidance and DME, as well as the signal processing for the 
angle data. The accuracy of the angle guidance information is preserved 
even under heavy multipath conditions by the use of a digitally implemented, 
matched tracking filter that acquires and tracks the direct signal. In 
acquiring the angle data, the processor employs a search algorithm that 
prevents lock-on to bright flashes or other spurious signals. Once in 
track, the receiver verification circuitry continuously checks the video 
spectrum to assure that the tracked signal is the correct one. In case of 
failure to verify, the receiver is forced to re-acquire the signal. 
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Table F-l. Functional Characteristics of the FAA-K Airborne Equipment 


Item 


Description 

Guidance Functions 

Azimuth, Elevation 1, Elevation 2, Back 


Azimuth, DME, 

Aux Data 

Coverage 

AZ, EL, DME 

BAZ EL2 

Horizontal 

±60 deg 

±40 deg ±45 deg 

Vertical 

1 to 

20 deg 1 to 8 deg 

Range 

30 nmi 

5 nmi 5 nmi 

MGA 

50 feet or TD 


Accuracies (2 CT) 

MGA 

Wide Angle or 
Max Range 

Noise: Azimuth 

0.026 deg 

0.10 deg 

Elevation 1 

0.06 deg 

0.08 deg 

Elevation 2 

0.03 deg 

- 

Back AZ 

0.052 deg 

0.052 deg 

DME (l<r) 

±20 ft 

±20 ft ±0.1%R 

Bias: Azimuth 

0.03 deg 

0.06 deg 

Elevation 1 

0.045 deg 

0.06 deg 

Elevation 2 

0.025 deg 

- 

Back AZ 

0.06 deg 

0.12 deg 

DME <1 <T) 

+20 ft 

±20 ft ±0 . 17.R 

Airborne Antenna Coverage 



Horizontal 


360 deg 

Vertical 

+5, 

-40 deg 


AZ ELI 

EL2 BAZ DME 

Acquisition Time (sec) 

1.6 1.4 

1.4 1.6 .25 

Verification Interval 

continuous 

(1 sec delay) N/A 
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Table F-2. Specific Performance Features of Configuration K, FAA 
Airborne Equipment 


Item 

Description 

Lateral Path 

Azimuth Select 

Pilot select of ±30 deg in 5-degree increments 

Sensitivity 

ILS compatible deviation with course softening 
option within desired range. 

Wide Angle 

±60 degree suitable for display. 

Missed Approach 

Automatic front-to-back AZ switching with 
DME. 

Vertical Path 


Glidepath Select 

2 to 12 degrees in 0.5-degree increments 

Sensitivity 

ILS compatible, with course softening option 
within desired range. 

Coordinates 

Conical, with option of ELI planar equivalent 
through DME algorithm. 

Flare Altitude and Rate 

Suitable for Collins 860F-1 display and 
11SA435 A flare coupler. 

Range and Range Rate 

Suitable for Collins 860-3 digital DME 
indicator. 

Monitoring 

On-line monitoring, self-monitoring push-to- 
test confidence test. Comparison test. 

Altitude/Range Discretes 

Marker beacon, flare, and decrab options. 
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The integrity of the angle receiver is enhanced by a combination of on-line 
monitoring, self-monitoring, and a press-to-test confidence check. The on- 
line monitoring is achieved by requiring the correct function identity (FI) 
to be decoded each data frame, and by the requirement that a minimum signal 
level be present before decoding of FI can start. The self-monitoring design 
makes extensive use of the microprocessor to monitor and test the status of 
the unit in a time share mode to ensure that the MLS is operational. The 
press-to-test confidence check circuitry injects a video signal to produce 
a cross-pointer deviation of predetermined magnitude and direction which 
provides a check, not only of the angle processor, but also of the interface 
circuits and pilot display. 

The HR-800 DME interrogator shares the channel selection local oscillator and 
RF front end with the angle data receiver. The channel selection also provides 
the pulse pair coding of the DME transmitter and decoding of the DME 
transponder signals. The DME range information may be used in automatic course 
softening of angle deviation signals, as well as in the computation of altitude 
to provide selectable Decision Height (DH) annunciation. 

The HN-400 control unit panel provides keyboard entry selection of any of 
200 MLS angle channels and the paired DME channel. The panel also allows key- 
board selection of azimuth elevation paths for deviation output reference. 

The selected paths, as well as the selected channel, are full-time displayed 
by gas-discharge type characters. The panel also has the confidence test 
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pushbutton and MLS subsystem status lights which are activated by the 
self-monitoring circuitry in the MLS equipment. The HN-400 control panel 
further provides access to the HN-700 microprocessor programming, wherein 
certain operatim al options may be implemented. For example, although 
normal Elevation 1 output is conical, both deviation and total angle may 
be converted to a planar equivalent through keyboard entry of proper 
coding to implement an appropriate DME conversion algorithm. 

The HL-362 and HL-363 omnidirectional antennas are provided to assure 
full MLS airborne antenna coverage during the diverse aircraft maneuvers 
associated with curved approach paths, missed approach, and departures. 
The HL-181 sector horn antenna provides effective gain enhancement of the 
guidance signals, as well as assuring adequate coverage for the critical 
final approach phase of the terminal area mission. 
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